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Chapter  1 
Introduction 


This  volume  gives  a  summapy  the  major  technical  achievements  on  Con¬ 
tract  No.  F33615-86-K-1023  for  the  period  April  1986  to  December  1989. 
The  overall  goals  of  this  contract  are  to  study  and  determine  solutions  for 
bistatic  first  order  and  higher  order  Uniform  Geometrical  Theory  of  Diffrac¬ 
tion  terms,  that  provide  an  accurate  and  efficient  means  to  calculate  high 
frequency  scattering  from  large  complex  geometries.  Algorithms  have  been 
developed  to  incorporate  these  techniques  into  a  user  oriented  computer 
code  referred  to  as  the  RCS  -  Basic  Scattering  Code  (RCS-BSC)  (1]. 

This  volume  presents  discussions  of  theoretical  studies  that  have 
been  completed  during  this  time  frame.  Chapter  2  discusses  various  meth¬ 
ods  for  determining  the  first  order  scattering  from  flat  plate  structures. 
It  presents  a  newly  developed  far  zone  comer  (vertex)  diffraction  coeffi¬ 
cient.  Complete  details  are  given  in  References  [2,3].  Chapter  3  presents 
a  discussion  of  a  form  of  edge  wave  -  vertex  interaction.  It  is  specifically 
formulated  in  this  case  for  a  source  excitation.  This  is  an  intermediate  step 
for  determining  the  far  zone  edge  wave  solution.  More  details  are  given 
in  References  [4,5].  The  extension  of  this  work  to  the  far  zone  edge  wave 
solution  is  presented  in  Chapter  4.  Other  higher  order  double  diffraction 
interactions  across  the  flat  plate  faces  are  discussed  in  Chapter  5. 


1 


Chapter  2 

Comparison  of  Methods  for 
Far  Zone  Scattering  from  a 
Flat  Plate  and  Cube 


2.1  Introduction 

The  validity  of  vaiious  methods  for  determining  the  far  zone  bistatic  scat¬ 
tering  from  a  flat  plate  and  convex  flat  plate  structure  such  as  a  cube  is 
presented  in  this  chapter.  This  is  accomplished  by  comparing  the  meth¬ 
ods  in  various  bstfic  situations.  The  specific  techniques  to  be  compared 
in  this  study  are  the  classical  equivalent  currents  with  “stripping”  [6],  the 
previous  corner  diffraction  coefficient  [6],  the  newly  developed  equivalent 
currents  by  Michaeli  [7],  and  an  extension  to  this  method  cast  in  the  form  of 
a  Uniform  Geometrical  Theory  of  Diffraction  (UTD)  far  zone  corner  diffrac¬ 
tion  coefficient  [2].  In  addition,  the  Method  of  Moments  (MOM)  using  the 
Electromagnetic  Surface  Patch  (ESP)  code  [8]  and  measurements  from  The 
Ohio  State  University  ElectroScience  Laboratories  compact  range  are  used 
to  further  validate  the  results. 

A  recent  paper  by  Ludwig  [9]  compares  three  methods  for  backscattering 
from  a  cube,  that  is,  the  MOM  using  the  Numerical  Electromagnetics  Code 
(NEC-MOM),  physical  optics  (PO),  and  the  previous  UTD  comer  diffrac¬ 
tion  solution.  In  this  chapter,  it  will  be  shown  that  methods  which  give 
comparable  results  for  backscatter  can  differ  for  bistatic  scattering.  The 
emphasis  here  is  to  present  basic  examples  that  can  be  used  to  validate 
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existing  codes  and  to  suggest  a  numerically  efficient  and  accurate  method 
for  convex  flat  plate  structures  to  first  order. 

An  approximate  expression  for  the  far  zone  field  scattered  by  the  vertex 
of  a  finite  perfectly  conducting  wedge  is  presented  in  this  regard.  The  solu¬ 
tion  is  cast  in  the  form  of  the  UTD  and  is  based  on  asymptotic  equivalent 
currents  found  using  modified  PTD  concepts  (7,2).  The  faces  of  the  wedge 
must  be  flat  (the  normal  to  each  individual  face  is  a  constant  everywhere  on 
the  face  except  at  the  edge)  and  the  edges  must  be  straight.  For  plane  wave 
incidence  from  an  arbitrary  direction,  the  first  order  contribution  from  each 
vertex  to  the  far  zone  scattered  field  is  obtiuned. 

Since  diffraction  is  a  local  phenomena  at  high  frequencies  the  results 
obtuned  for  a  finite  wedge  may  be  applied  to  much  more  complex  bodies 
made  up  of  simple  shapes.  The  field  scattered  by  a  three-dimensional  shape 
constructed  from  flat  plates  may  be  approximated  to  first  order  as  the  sum 
of  the  contributions  from  each  individual  corner.  The  first  order  solution 
should  be  reasonably  accurate  in  or  near  the  specular  regions  as  long  as 
the  object  is  convex.  A  convex  body  is  defined  here  as  a  closed  surface 
made  up  of  flat  plates  such  that  sdl  of  the  exterior  wedge  singles,  taken 
between  faces  and  exterior  to  the  surface,  are  greater  than  180  degrees.  A 
simple  example  of  an  object  that  does  not  meet  this  requirement  is  a  comer 
reflector.  In  this  case,  the  effect  of  the  interaction  between  the  faces  must 
be  taken  into  account.  Higher  order  effects  such  as  double  diffraction  [10] 
and  edge  waves  [11]  are  not  considered  in  this  chapter. 

Note  that  the  restdts  presented  in  this  chapter  are  for  a  parallel  ray 
type  solution,  that  is,  for  a  radar  cross  section  result.  The  NEC  -  Basic 
Scattering  Code  (NEC-BSC)(12]  is  a  near  zone  formulated  code,  that  is  it 
has  a  finite  range  involved.  The  UTD  solutions  are  slightly  different  for  this 
non-parallel  ray  case.  The  capabilities  of  the  NEC-BSC  and  a  comparable 
far  zone  code  called  the  RCS-BSC  are  discussed  in  Reference  [13]. 


2.2  Theoretical  Background 

There  are  many  approximate  solutions  to  the  scattered  field  from  a  finite 
perfectly  conducting  wedge.  Physical  Optics  Bind  its  extension  the  Physical 
Theory  of  Diffraction  [14]  is  surface  and  edge  current  based.  Geometri¬ 
cal  Optics  (GO)  and  its  extensions  the  Geometrical  Theory  of  Diffraction 
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(GTD)  [15]  and  the  Uniform  Geometrical  Theory  of  Diffraction  [16]  are  ray 
based.  The  Method  Equivalent  Currents  (MEG)  [17]  is  an  intermediate 
type  solution  that  was  developed  to  handle  caustic  regions  in  the  GTD. 
This  has  been  augmented  with  the  concept  of  stripping  to  provide  better 
answers  for  flat  plate  problems  [6].  Recently,  Michaeli  [18]  showed  a  mere 
rigorous  approach  in  deriving  equivalent  currents.  This  was  shown  to  be  re¬ 
lated  to  the  incremental  length  method  of  Mitzner  [19]  by  Knott  [20].  These 
equivalent  currents  still  had  singularity  problems  that  have  been  remedied 
by  Michaeli  [7]  using  a  skewed  coordinate  system.  Uflmtsev  also  derived  a 
similar  solution  [21,22]. 

The  above  solutions  can  be  cast  in  a  comer  diffraction  coefficient  form. 
These  UTD  ray  type  solutions  have  the  advantage  of  being  efficient  for 
far  zone  flat  plate  problems  since  only  the  fields  scattered  from  the  comers 
need  to  be  added.  It  also  has  the  advantage  that  the  results  correlate  to  the 
scattering  centers  seen  in  high  resolution  measurements.  Just  the  comer 
diffraction  coefficient  forms  are  outlined  in  this  section. 

A  previous  diffraction  coefficient  for  a  comer  formed  by  the  intersection 
of  two  straight  edges  was  derived  by  Burnside  and  Pathak  [6].  It  is  based  on 
the  asymptotic  evaluation  of  the  radiation  integral  containing  the  equivalent 
currents  of  Ryan  and  Peters  [17].  The  result  was  then  empirically  modified 
so  that  the  diffraction  coefficient  would  not  change  sign  abruptly  as  it  passes 
through  the  false  shadow  boundaries.  It  was  derived  for  spherical  wave 
incidence  and  remains  valid  for  cases  when  the  diffraction  point  is  near  the 
corner  since  the  integral  was  evaluated  for  a  saddle  point  near  an  end  point; 
however,  only  the  far  zone  result  is  shown  here.  The  comer  diffracted  field 
due  to  one  comer  and  one  edge  in  the  case  of  plane  wave  incidence  and  a 
far  zone  receiver  is  ^ven  by 
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Figure  2.1:  Definition  of  angles  for  the  Previous  Comer  Diffraction  Coeffi¬ 
cients. 
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where  is  the  integer  which  most  nearly  satisfies  2nnN^  —  =  T’t*  *nd 
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where  the  angles  are  shown  in  Figure  2.1.  The  sign  on  the  diffraction  co¬ 
efficient  may  be  plus  or  minus  depending  on  which  endpoint  of  the  edge 
is  being  considered.  The  correct  sign  in  front  of  the  (7,,h  terms  in  Equa¬ 
tion  (2.2)  is  chosen  based  on  the  direction  of  the  edge  vector  shown  in 
Figure  2.1. 

It  is  assumed  that  the  incident  field,  and  therefore  the  scattered  field, 
is  a  time  harmonic  field  with  time  dependence  given  by  which  is  sup¬ 
pressed. 
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The  new  far  zone  corner  diffraction  solution  is  based  on  the  PTD  and 
cast  into  the  form  of  the  MEC  and  then  into  a  UTD  diffraction  coefficient. 
This  is  done  as  follows.  The  PO  is  first  used  to  approximate  the  currents 
resulting  in  a  double  integral  over  the  surface.  Stokes  theorem  is  then  ap¬ 
plied  to  reduce  the  equation  to  a  line  integral  [23,24].  The  Michaeli  currents 
are  added  to  produce  a  total  first  order  MEC  resiilt.  This  integral  is  then 
evaluated  using  the  method  of  stationsiry  phase  to  obtain  the  contribution 
from  each  comer  [25]. 

The  new  corner  diffraction  coefRcients  are  given  in  a  form  similar  to 
previous  expressions  for  diffraction  coefficients: 
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where  the  plus  or  minus  sign  is  chosen  depending  on  which  endpoint  con¬ 
tribution  is  being  calculated.  The  minus  sign  is  used  for  the  corner  con¬ 
tribution  associated  with  the  negative  t-axis,  while  the  plus  sign  is  used 
for  the  corner  contribution  associated  mtb  the  positive  t-axis.  The  edge 
fixed  coordinates  shown  in  Figure  2.2  are  chosen  such  that  n  is  the  outward 
normal  of  the  0-face,  t  is  tangent  to  the  edge,  the  positive  b-axis  lies  on 
the  0-face,  and  t  =  6  x  n.  The  expressions  for  and 
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where  the  +  sign  is  associated  with  ■>  *nd 

while  the  —  sign  is  associated  with  the  and  df*^  terms.  The 

angles  are  defined  in  Figure  2.2.  Since  only  convex  structures  are  consid¬ 
ered  here  proper  shadowing  of  the  rays  is  fairly  simple.  The  shadowing  of 
the  incident  field  is  accounted  for  by  E^,  and  E^,  which  are  the  components 
of  the  GO  incident  field.  The  shadowing  of  the  diffracted  ray  is  more  com¬ 
plicated.  The  contributions  from  the  LPO  and  PO  components,  dj;^^  and 
are  present  everywhere.  The  UTD  components,  djfj^,  are  shadowed 
like  diffracted  fields.  They  do  not  contribute  if  the  observation  point  is 
inside  the  wedge  {<f>  >  nir). 

For  the  special  case  of  a  fiat  plate  (n  =  2)  the  contribution  from  both 
faces  may  be  found  using 
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Figure  2.2:  Definition  of  the  Angles  used  in  the  New  Corner  Diffraction 
Coefficients. 
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where  7,  a,  and  the  other  variables  have  been  defined  previously. 

It  is  interesting  to  note  that  by  writing  the  equations  for  the  Michaeli 
equivalent  ctirrents  and  the  new  comer  diffraction  coefficients  in  cotangent 
form  provides  more  insight  into  the  connection  of  the  new  solutions  with 
the  previous  methods.  The  new  parameters  separate  out  the  optics  currents 
and  diffraction  currents.  This  separation  manifests  itself  in  new  parameters 
for  the  <f>  angles.  They  arise  from  the  asymptotic  evaluation  of  the  currents 
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Figure  2.3:  Two  wavelength  plate  in  the  x-z  plate. 

in  the  skewed  coordinate  system  chosen  in  physically  meaningful  directions. 
The  LPO  factor  (7)  is  related  to  the  projection  of  the  average  of  the  incident 
and  diffraction  planes  on  to  the  plane  of  the  plate  [2].  The  PO  and  UTD 
factor  (a)  relates  to  the  projection  of  the  Keller  diffraction  cone  on  to  the 
plane  of  the  plate.  It  is  easy  to  see  in  this  form  that  in  the  Keller  directions 
the  LPO  and  PO  cancel,  leaving  the  UTD  result  formally  used  in  many 
solutions. 


2.3  Comparisons 

The  first  example  compares  the  Ryan  and  Peters  equivalent  currents,  the 
previous  and  new  corner  diffraction  solutions.  The  simple  example  of 
backscatter  from  a  two  wavelength  square  plate  lying  in  the  x-z  plane, 
as  shown  in  Figure  2.3,  is  used.  This  illustrates  that  for  backscatter  these 
different  methods  produce  very  similar  results,  except  for  the  very  low  level 
regions. 

The  co-polarized  fields,  in  the  principal  plane,  calculated  using  the  three 
different  methods  are  shown  in  Figures  2.4  and  2.5.  All  three  methods 
give  essentially  the  same  results  for  the  principal  plane  pattern  cuts  shown 
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Figure  2.4:  Backscatter  from  a  2  wavelength  plate  (9  =  90*  pattern). 

here.  This  is  not  surprising  since  the  major  contributions  to  the  fields 
are  the  scattering  from  the  two  edges  in  their  KeUer  cone  directions.  The 
new  comer  diffraction  solution  reduces  to  the  Ryan  and  Peters  equivalent 
current  solution  for  paints  on  the  Keller  cone  [2],  and  the  previous  comer 
diffraction  solution  is  essentially  the  same  as  Ryan  and  Peters  equivalent 
current  solution  for  most  regions  of  space.  The  results  in  Figure  2.5  are  for 
the  horizontal  (<r^)  polarization.  For  a  knife  edged  plate  such  as  this,  the 
scattered  field  should  be  zero  in  the  plane  of  the  plate.  Note  that  this  is 
not  the  case  in  these  first  order  results.  The  higher  order  terms  (i.e.  the 
double,  triple  etc.  diffractions)  produce  the  null  for  this  polarization  when 
they  are  included. 

For  patterns  away  from  the  principal  plane,  the  higher  levels  are  the 
same  but  the  lower  levels  differ.  This  is  illustrated  by  taking  a  conical  cut 
(B  =  60”)  for  the  two  wavelength  plate.  The  results  for  the  same  three 
methods  used  previously  are  shown  in  Figures  2.6  and  2.7.  In  this  case 
the  methods  agree  well  for  the  main  lobe,  however,  they  differ  in  the  lower 
levels  of  the  pattern. 

The  differences  in  the  three  methods  mentioned  earlier  are  greatly  in¬ 
creased  for  bistatic  scattering  problems.  The  bistatic  scattering  from  a 
square  plate  two  wavelengths  on  a  side  is  examined  to  illustrate  this  point. 
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Figure  2.5:  Backscatter  from  2  wavelength  plate  (9  =  90“  pattern). 
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Figure  2.6:  Backscatter  from  2  wavelength  plate  (9  =  60“  pattern). 
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Figure  2.7:  Backscatter  &om  a  2  wavelength  plate  {6  =  60°  pattern). 

The  complete  scattering  matrix  (all  four  values  of  <r)  is  found  for  s  plate  in 
the  x-y  plane  with  a  fixed  source  located  at  =  45°  and  ^‘  =  0°  as  shown 
in  Figure  2.8.  The  results  for  the  ^  =  60°  pattern  cut  are  compared  with 
the  previous  comer  diffraction  solution  and  Method  of  Moment  calculations 
for  co-polarized  fields  in  Figures  2.9  and  2.10.  Similarly  the  results  for  the 
cross-polarized  fields  are  given  in  Figure  2.11  and  Figure  2.12.  Overall  the 
new  solution  agrees  well  with  the  Method  of  Moment  calculations  and  does 
not  exhibit  the  discontinviities  which  appear  near  0  =  240°  and  0  =  300° 
in  the  previous  comer  diffraction  solution.  The  discontinuities  in  the  pre¬ 
vious  corner  diffraction  solution  are  caused  by  the  so  called  false  shadow 
boundaries  where  the  associated  two-dimensional  problem  passes  through 
a  shadow  boundary,  but  the  three-dimensional  problem  in  reality  does  not. 
The  Ryan  and  Peters  equivalent  current  results  are  not  shown  here,  but 
they  behave  differently  for  rimilar  reasons;  that  is,  the  solution  still  con¬ 
tains  two  dimension  information  in  regions  that  it  should  not.  In  the  region 
from  0  fs  60°  to  120°  (i.e.  near  the  plane  of  the  plate)  the  new  solution  and 
the  Method  of  Moments  solution  differ  by  more  than  10  dB.  It  is  suspected 
that  most  of  these  differences  are  due  to  the  effects  of  higher  order  terms 
(double  and  triple  diffraction,  edge  waves)  which  are  not  included  in  the 
new  solution. 
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Figure  2.8:  2A  square  plate  in  the  x-y  plane  with  a  fixed  source  at  S'  =  45“ 
and  =  0. 


Figure  2.9:  Co-polarized  RCS  in  the  ^  =  60“  plane  of  a  2A  square  plate 
with  a  ^  polarized  fixed  source  at  6*  =  45“,  ^’  =  0“. 
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Figure  2.10:  Co-polarized  RGS  in  the  (f>  —  60°  plane  of  a  2A  square  plate 
with  a  polarized  fixed  source  at  =  45°,  =  0®. 


Figure  2.11:  Cross-polarized  RCS  in  the  ^  =  60®  plane  of  a  2A  square  plate 
with  a  polarized  fixed  source  at  ^  =  45®,  ^‘  =  0®. 
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Figure  2.12:  Cross-polarized  RCS  in  the  ^  =  60°  plane  of  a  2A  square  plate 
with  a  polarized  fixed  source  at  B*  =  45°,  =  0°. 

In  this  example  the  new  solution  is  compared  to  backscatter  measure¬ 
ments  [26]  made  at  10  GHz  on  a  6”  cube.  The  geometry  of  the  cube, 
tilted  45°  in  the  x-z  plane,  is  illustrated  in  Figure  2.13.  The  results  for 
the  H-plane  and  E-plane  patterns  taken  in  the  x-y  plane  are  given  in  Fig¬ 
ures  2.14  and  2.15,  respectively.  The  results  agree  well  to  first  order  over 
most  regions  of  the  pattern.  The  discrepancies  are  probably  due  to  a  com¬ 
bination  of  higher  order  terms  not  being  included  in  the  analysis  and  in 
measurements  errors.  The  error  in  the  measurements  is  likely  two  fold. 
First  the  faces  of  the  cube  model  were  misaligned  slightly  so  they  did  not 
form  edges  as  sharp  as  may  be  required.  Secondly,  it  seems  that  there  was 
some  deviation  from  the  desired  pattern  cuts  as  can  be  seen  from  the  lack  of 
symmetry  in  the  measured  patterns.  In  any  case,  they  confirm  the  validity 
of  the  new  comer  diffraction  solution  within  first  order  accuracy  for  wedge 
type  structures. 


2.4  Discussion 

The  new  corner  diffraction  coefficient  in  the  above  examples  has  been  shown 
to  provide  improved  results  over  other  methods,  especiaJly  in  bistatic  Bit- 


15 


Figure  2.16:  RCS  for  the  B  =  89“  cut  of  a  2A  square  plate  with  a  polarized 
fixed  source  at  B*  =  45®,  =  0®. 


uations.  The  Michaeli  equivalent  currents  have  not  been  shown  since  they 
provide  the  same  results  as  the  new  comer  diffraction  coefficient.  Certain 
properties  of  these  new  solutions,  however,  may  stiU  cause  patterns  taken 
in  some  repons  of  space  to  be  discontinuous. 

It  has  been  shown  [7,2]  that  £)|  and  Df^  do  not  tend  to  definite  limits 
as  i  — »  O'  (i.e.  the  intersection  of  the  associated  half-plane  and  the  Keller 
cone),  where  &  =  tain 0*  +  6 cos/?',  but  they  remain  bounded.  In  practice, 
this  means  that  both  Z),  and  and  therefore  Ep  and  E^,  are  discontin¬ 
uous  at  this  point  in  the  pattern.  A  simple  example  illustrates  how  this 
discontinuity  can  affect  a  pattern.  The  bistatic  RCS  from  the  flat  plate 
shown  earlier  in  Figure  2.8  is  considered.  The  source,  linearly  polarized 
in  the  direction,  remains  fixed  at  (?*  =  45®  and  =  0®  while  the  pat¬ 
tern  is  taken  near  the  x-y  plane  (B  =  89®).  The  bistatic  RCS  is  given  in 
Figures  2.16  and  2.17  for  the  co-polarized  and  cross-polarized  fields,  re¬ 
spectively.  The  abrupt  null  at  ^  w  135®  in  the  co-polarized  pattern  and 
the  spike  at  the  same  location  in  the  cross-polarized  pattern  are  due  to 
discontinuities  in  the  contribution  from  edge  4  (indicated  in  Figure  2.8). 
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Figure  2.17:  RCS  for  the  B  —  89°  cut  of  a  2A  squau'e  plate  with  a  polarized 
fixed  source  at  ^  =  45°,  =  0°. 


The  point  <l>  w  135°  coincides  with  (3^  =  0'^  and  04  ss  0  where  0^,  13'^,  and  04 
are  the  edge  fixed  coordinates  for  edge  4.  Due  to  the  geometry  0^  fa  ^  and 
04  %  d  so  the  discontinuity  in  is  due  to  the  discontinuity  in  and, 
likewise,  the  discontinuity  in  is  due  to  the  discontinuity  in  D^. 

Therefore,  the  discontinuity  in  the  new  diffraction  coefficients  at  the 
intersection  of  the  Keller  cone  and  the  infinite  half  plane  associated  with 
the  edge  =  0'  and  0  =  0)  may  be  expected  to  cause  discontinuities 
or  narrow  spikes  depending  on  the  polarization  and  the  pattern  cut.  As 
the  examples  illustrate  these  disturbances  only  affect  a  typical  pattern  cut 
for  around  5°  to  10°.  In  addition,  they  are  in  the  low  level  regions  of  the 
returns. 

It  is  easily  seen  that  the  diffraction  coefficients  and  are  discon¬ 
tinuous  as  the  source  passes  through  the  half  plane  0'  =  tt.  In  the  general 
case  of  bistatic  scattering,  these  discontinuities  in  the  sign  of  the  field  scat¬ 
tered  by  a  corner  will  result  in  discontinuities  in  the  total  scattered  field. 
However,  the  diffraction  coefficients  are  continuous  here  (0'  =  ir)  for  the 
special  case  of  backscatter. 
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2.5  Conclusions 


The  objective  of  this  chapter  has  been  to  compare  different  methods  for 
the  analysis  of  the  high  frequency  far  zone  scattering  from  flat  plate  and 
convex  flat  plate  type  structures.  Ryan  and  Peters  equivalent  currents  and 
the  previous  corner  diffraction  coefficient  are  compared  with  the  Michaeli 
equivalent  currents  and  the  new  corner  diffraction  coefficient.  The  method 
of  moments  and  measurements  are  also  used  to  validate  the  solutions. 

It  has  been  shown  that  for  backscatter  all  the  methods  compare  rea¬ 
sonably  within  engineering  accuracy.  For  bistatic  scattering,  however,  the 
two  dimensional  nature  of  the  old  methods  lead  to  inaccuracies.  The  newer 
methods,  based  on  more  rigorous  three  dimensional  analysis,  remove  most 
of  these  problems. 

A  new  comer  diffraction  coefficient  is  presented  that  provides  an  efficient 
and  accurate  scdution  to  within  first  order.  It  provides  the  same  level  of 
accuracy  as  the  Michaeli  equivalent  currents  with  the  added  benefit  of  not 
needing  integrations  for  flat  plates.  All  the  optics  and  edge  scattering 
effects  have  been  lumped  into  the  comers  of  the  plate  with  nice  physical 
interpretations. 
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Chapter  3 


Diffraction  of  Dipole  Excited 
Edge  Waves 

3.1  Introduction 

The  diffraction  of  acoustical  waves  by  a  conducting  elliptic  cone  and  a 
quarter  plane  were  studied  by  Krauss  and  Levine  [27]  and  Radlow  [28,29], 
respectively.  Satterwhite  and  Kouyoumjian  [30]  analyzed  the  vector  elec¬ 
tromagnetic  problem  and  developed  a  Green’s  dyadic  for  a  source  radiating 
in  the  presence  of  an  angular  sector.  However,  their  solution,  expressed 
in  terms  of  non-closed  form  Lame  functions,  is  cumbersome  for  numerical 
calculations.  Furthermore,  it  has  not  so  far  appeared  possible  to  asymp¬ 
totically  identify  a  “tip  diffraction  coefficient”  from  this  eigenfunction  rep¬ 
resentation. 

Burnside  and  Pathak  proposed  in  [31]  a  corner  diffraction  coefficient 
which  successfully  predicted  the  corner  effect  of  several  plate  structures. 
Their  solution  is  based  on  the  asymptotic  evaluation  of  the  radiation  in¬ 
tegral  involving  the  equivalent  currents  that  would  exist  in  the  absence  of 
the  comer.  A  corner  diffraction  term  is  then  established  by  empirically 
modifying  the  final  result.  Sikta  [11]  applied  a  limiting  process  to  identify 
the  wave  diffracted  by  the  corner  and  propagating  along  one  of  the  edges 
of  a  polygonal  plate.  Further,  introducing  an  empirically  established  “re¬ 
flection  coefficient”,  he  utilized  his  expression  to  the  calculation  of  multiple 
diffractions  between  adjacent  corners  of  a  fiat  plate. 

In  this  Chapter,  the  vertex  and  edge  diffraction  of  an  electromagnetic 
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wave  gtiided  along  one  of  the  edges  of  a  semi-infinite  wedge  is  studied 
asymptotically.  Explicit  expressions  of  the  edge  wave  are  presented  in  Sec¬ 
tion  3.2,  deduced  from  the  limiting  behavior  of  the  Green’s  dyadic  for  an 
infinite  wedge  [32,33]. 

A  first  approximation  of  the  vertex  diffracted  field  is  developed  in  Sec¬ 
tion  3.4,  based  on  the  radiation  of  the  surface  current  that  would  be  induced 
by  the  edge  wave  for  an  infinite  wedge,  which  is  subsequently  truncated. 
The  asymptotic  evaluation  of  the  surface  current’s  radiation  integral  ap¬ 
propriately  encoimters  the  edge  wave  singularity,  consistent  with  Meixner’s 
edge  condition.  To  the  field  expressions  thus  derived,  a  fringe  current  effect, 
which  is  asymptotically  incorporated  in  the  radiation  integral  of  Michaeli’s 
fringe  edge  equivalent  currents  [7],  is  superimposed.  Edge  wave  vertex  and 
edge  diffraction  coefficients  can  then  be  established  from  the  asymptotic 
field  expressions.  The  pertinent  analyms  is  developed  in  Section  3.5. 

The  validity  of  the  approach  is  confirmed  via  comparison  with  Method 
of  Moments  results  and  pattern  measurements  for  a  small  dipole  radiating 
in  the  close  vicinity  of  one  of  the  edges  of  a  polygonal  plate. 


3.2  Edge  Waves 

The  term  **edge  waves”  in  the  present  work  defines  waves  stron^y  coupled 
with  the  edge  of  a  wedge.  The  edge  wave  is  a  Maxwellian  field  guided  by 
the  edge,  exhibits  the  proper  edge  singularity  and  can  be  excited  either 
by  a  plane  wave  at  grazing  incidence,  or  by  a  dipole  radiating  in  the  close 
vicinity  of  the  edge.  The  vertex  of  a  polyhedron  illuminated  by  a  plane  or 
spherical  field  can  also  excite  an  edge  wave.  Independent  of  the  excitation, 
but  sufficiently  far  from  its  source,  the  edge  wave  is  and  can  be  treated 
as  a  ray  optical  field,  i.e.,  its  phase  is  proportional  to  the  distance  from 
the  source  (reference)  to  the  observation  point  and  its  amplitude  varies  in 
such  a  way  that  the  energy  within  a  tube  originated  at  the  source  point  is 
preserved.  Further,  its  Poynting  vector  is  radially  directed.  Application  of 
ray  optical  techniques  is,  however,  not  straightforward,  mainly  due  to  the 
singular  behavior  of  the  paraxial  fields. 

Let  us  consider  the  canonical  problem  of  a  point  electric  source  j*  = 
p*  S(f  —  r')  that  radiates  in  the  presence  of  a  perfectly  conducting  infinite 
wedge.  The  field  produced  by  this  dipole-wedge  configuration  is  formally 


22 


written 


E(f)  =  jufofiof(f-,  P) .  p%  (3.1) 

where  r(r;r')  is  the  Green’s  dyadic  for  the  wedge  [32].  For  convenience, 


Figure  3.1:  Dipole  excited  edge  wave. 

the  source  point  is  located  at  r'  =  p',  =  ir/2  in  the  reference  coordinates, 

with  the  z*axiB  being  coincident  with  the  edge  and  the  s-axis  parallel  to 
face  1  of  the  wedge  (Fig.  3.1).  The  related  system  of  spherical  coordinates 
is  (r,  /?o,  4i).  The  limiting  form  of  the  field  in  eq.  (3.1)  for  small  values  of  the 
parameter  e  =  kp' sin  0o  is  of  interest  here.  Specifically,  the  Green’s  dyadic 
assumes  the  power  series  expansion 

=  e-'  f;  (3.2) 

m=0  n=0 

where  1/  =  ir/7  and  7  is  the  external  wedge  angle.  The  dyadic  terms 
^rnni^i  p' ,  ^  =  0»  1»  •  •  •  i  1>  •  •  •  >  b®  derived  in  a  straightfor¬ 

ward  manner,  if  one  substitutes  the  series  expansion  of  the  Hankel  functions 
in  the  eigenfunction  representation  of  r(r;r').  At  the  limit  e  — »  0,  one  can 
retain  the  leading  term  in  eq.  (3.2)  to  approximate  the  total  field  which, 
for  ifcr  >  >  1 ,  is  given  by 
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where 


and 


(3.3) 

(3.4) 

(3.5) 


is  a  constant  source  factor.  Eq.  (3.3)  defines  an  edge  guided  wave,  which 
dominates  in  the  paraxial  region,  revealing  the  strong  coupling  between  the 
dipole  source  and  the  edge.  It  is  a  spherical  wave  that  satisfies  Maxwell’s 
equations,  the  boundary  conditions  on  the  perfectly  conducting  surface,  as 
well  as  Meixner’s  edge  condition. 

In  the  case  of  a  plane  wave  at  grazing  incidence,  which  is  essentially 
a  quasi-two  dimensional  situation,  one  can  implement  a  similsu'  limiting 
procedure  on  the  cylindrical  harmonic  expansion  of  the  total  field  [33]  to 
extract  the  edge  wave.  In  particular,  with  —*  0,  the  dominant  term  in 
the  series  expansion  of  the  total  field  in  the  nrighborhood  of  the  edge  equals 


=  -AirC{v)  4>)ikpY~^  {p  «n  ^  cos  u4>)  txp{-jkz),  (3.6) 


where  S{0',(l>')  =  sin"  ^  0\Eop'  amu<f)'  +  Eo^  cosi/^')  is  again  an  excitation 
related  factor. 


Both  forms  (3.3)  and  (3.6),  although  derived  from  n  limiting  rather  than 
asymptotic  estimation,  agree  in  the  paraxial  region  with  the  leading  term 
in  the  asymptotic  analyris  of  Pearson  [34]. 


3.3  Diffraction  formulation 

The  study  of  edge  guided  waves  at  an  infinite  wedge,  outlined  in  the  previ¬ 
ous  section,  is  not  adequate  to  describe  the  radiation  of  a  dipole  in  the  close 
vicinity  of  one  of  the  edges  of  a  realistic  fiat  plate  structure.  For  instance, 
regarding  the  simple  trihedral  configuration  depicted  in  Fig.  3.2,  in  addi¬ 
tion  to  a  direct  field  associated  with  the  interaction  of  the  dipole  and  the 
neighboring  edge,  diffracted  rays  initiated  at  the  tip  and  the  adjacent  edges 
should  be  introduced  in  an  accurate  high  frequency  solution.  Namely,  an 
improved  approximation  of  the  total  field  can  be  expressed  as 

+  ^  +  +  (3.7) 


24 


In  the  above  equation,  •  7  constitutes  the  dipole  and  guiding  edge 
interaction;  in  essence,  it  is  a  multiplicative  correction  of  the  edge  wave 
given  by  eq.  (3.3)  ,  by  the  transition  dyadic  7  to  be  determined,  thus 
assuring  a  finite  field  when  0^0,  i.e.,  at  the  paraxial  re^on  of  the  extension 
of  the  edge  in  free  space.  In  addition,  denote  edge  diffracted  fields, 
while  E‘  designates  the  ray  emanating  from  the  tip. 


FIELD  POINT 


In  principle,  if  j{r')  is  the  total  current  flowing  over  the  surface  S  of 
the  conducting  trihedron,  then  the  scattered  field  can  be  evaluated  via  the 
radiation  integral,  which,  in  the  Fresnel  or  FVaunhofer  region  of  the  surface, 
has  the  form 

E{s,)  =  jkZo  J  J^RxRx  j{P)  •  (3.8) 

Asymptotic  expansion  of  the  above  integral  by  the  method  of  stationary 
phase  [35],  in  terms  of  the  large  parameter  A;,  should  restilt  in  isolated 
mechanisms  associated  with  specific  parts  of  the  object  and,  in  our  case, 
corresponding  to  the  diffraction  terms  in  the  right  hand  side  of  eq.  (3.7). 
How.^ver,  use  of  the  radiation  integral  in  field  calculations,  requires  a  prior 
knowledge  (at  least  approximate)  of  the  current  j(r')  on  5. 
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Two  successive  estimates  of  the  current,  resulting  in  different  levels  of 
approximation  are  introduced.  First,  it  is  postulated  that,  on  the  faces  of 
the  trihedron  grazed  by  the  edge  wave,  j{r')  can  be  satisfactorily  approxi¬ 
mated  by  the  actual  induced  current  (edge  wave  current)  that  would  flow 
if  the  wedge  was  infinite,  while  it  is  presumed  zero  at  the  shadowed  face  of 
the  structure.  In  particular,  for  face  1  (Fig.  3.3(a)),  we  presume 

=  ft.  X  .  (3.9) 

where  is  the  edge  wave  at  an  infinite  wedge  and  ni  is  the  unit  vector 

normal  to  face  1. 

The  asymptotic  reduction  of  the  surface  integral  in  eq.  (3.8)  to  a  line 
integral  along  edges  (1)  and  (2)  in  Fig.  3.3(b)  reveals  an  equivalence  between 
the  radiation  of  the  edge  wave  surface  current  (eq.  (3.9))  and  the  Physical 
Optics  equivalent  currents  in  [7],  proportional  at  each  point  on  the  edge  to 
the  incident  edge  wave  (rather  than  a  plane  wave,  for  which  they  have  been 
originally  developed).  Based  on  this  observation,  to  the  truncation  effect 
which  results  from  eq.  (3.8),  a  non-uniform  (fringe)  current  effect  induced 
by  edges  (1)  and  (2)  is  superimposed  in  a  rather  indirect  way:  via  the 
radiation  of  Michaeli’s  fringe  equivalent  currents  for  edge  wave  incidence. 

In  a  consistent  rationale,  one  should  include  a  tip  excited  non-uniform 
current.  But  the  latter  can  not  be  extracted  from  existing  representations 
of  the  comer  diffraction  phenomenon  and  remains  unknown.  Numerical 
results  indicate,  however,  that  the  present  depiction  of  the  currents  is  suf¬ 
ficient  to  approximate  the  diffracted  edge  waves  at  the  vertex. 


3.4  The  radiation  of  the  edge  wave  current 

Let  us  consider  the  radiation  of  the  edge  wave  current  flowing  over  face 
1  of  the  trihedron.  The  results  for  face  2  can  be  derived  by  means  of  a 
simple  transformation.  One  may  introduce  the  oblique  system  of  coordi¬ 
nates  (u,<i)  associated  with  edge  (1)  in  Fig.  3.3(a)  and  defined  by  the  unit 
vectors 

u  =  i,  <1  =  X  sinai  —  i  cosoi, 

with  ai  denoting  the  angle  formed  by  the  guiding  edge  and  edge  (1).  Suf¬ 
ficiently  far  from  the  dipole  and  employing  the  results  of  Section  3.2,  one 
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Figure  3.3:  Geometry  of  the  edge  wave  diffraction  problem. 
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observes  that 


cxp( — 3 — o)  0  +  *  cos  tfo)» 

I’o 


(3.10) 


in  which 


To  =  y/ii  sin*  Oi  +  (u  -  ticosai  +  s')*, 

sin 00  =  ai/ro,  cos  d©  =  (u  —  tj  cos  Oj  +  s[.)fro, 

s'  is  the  distance  of  the  dipole  source  from  the  vertex  and  the  constant 
C{v)  is  defined  by  eq.  (3.4). 

The  vertex  diffracted  field  of  the  wedge  can  be  identified  asymptoti¬ 
cally  as  the  end  point  contribution  to  the  radiation  integral  in  eq.  (3.8). 
Classically,  employing  a  quadrature  expansion  of  the  phase  function,  one 
obtains 


£?r“'(-c) « - 


3k''Zo  sin^  ai  C(u)  A*(p*,  ^')  F{2kLc  sin*  f ) 


'  '  4irse(s0‘'+*(l -co8^) 

•  |(s' sin/9/^„(fc)  -  (sin^cosai  +  sinai  co8/3cos^)/2„_i(fc)]4 
-|-sinaisin^/^^_i(fc)^J .  (3.11) 

In  the  above  expression,  Hp{k)  signifies  the  end  point  contribution  to  the 
integral 

I-p{k)  =  f  exp  {-iA!(ii  +  ro)|«=o}  dt, 

Jo 

which  is  studied  in  Appendix  A. 

For  the  specific  semi-infinite  wedge  geometry  I?.p(k)  can  be  approxi¬ 
mated  by 


Ilf{k)  »  r(p)  k  P  exp(ipir/2)  exp{-ifc(sc  +  ^c)} 

_ i^p[fca?(*;;8;)] 

[nnoj  sin^ cos^  -I-  co8ai(l  —  cos/5)]*” 


in  which 


=  -•ign(»  -  Oi  -  A)  +«i) 


(3.12) 


(3.13) 
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and  are  the  distances  of  the  dipole  source  point  and  the  receiver  from 

the  origin  QJ  of  the  Keller  cone  of  diffracted  rays  initiated  at  edge  (1).  The 
branch  of  the  bracketed  expression  in  the  denominator  of  eq.  (3.12)  should 
be  chosen  according  to  (-1)'’  =  exp(jpn).  The  factor  F{-)  is  the  familiar 
edge  transition  function  of  the  Uniform  Theory  of  Diffraction  (UTD)  [16], 
while  the  edge  wave  transition  function  *  |)  defined  by 

FZp{x)  =  cxp{jpir/A)  (2*)*’/’  cxp{jzl2)  Z>_p[exp{>ir/4)v^],  (3.14) 

with  2?_p(-)  denoting  the  parabolic  cylinder  function  of  order  — p  (eq.  (A.4)). 
For  z  >  0  small  we  have 


Ffp(x)  «  v/5r  exp(jp7r/4)z*’^*  exp(jz/2) 

1  2exp(jV/4)  ^  i(p  - 1/2) 

[r(H±i)  r(*)  r(H±i) 

and  for  z  large 


(3.15) 


P(P+  1)  _  P{P  -H)(P  +  2)(p  +  3) 
Ajx  32z* 


(3.16) 


Retaining  only  the  dominant  term  in  the  expression  of  the  vertex  dif¬ 
fracted  wave  (dominant  with  respect  to  the  parameter  2)bs')  and  with  the 
approximation  [36] 


«i.a  +  «i.a  =«  («c  +  Jc)  -  ic  -!  ct's(cri,j 


.i-  a  M 


for  large  values  of  the  distance  parameter  =  «c«c/(*e+«c)»  obtains  the 
following  expression  for  the  vertex  diffracted  field,  in  terms  of  the  edges  (1) 
and  (2)  fixed  coordinate  systems  («„  ^\)  and  (Sc.  ^a).  respectively: 


^.ew/ -  X  ^  exp(jt/ar/2)  r(i/)  C(t/)  £^,(Q*^)  sin  u<l>'  -|-  cos  Vif)' 


sin 


■  cot  ^  F{2kLc  sin 


•"  2> 


sin’^tti  Ff,,[2feLeC08»  (g^)] 
(cOSQi  -I-  COB^iY 


sin*'  a,  FiJ^2kL^  cos’ 

(coSO]  -J-  COS/Jj)" 


I  exp(-jfese)  A 

J 


(3.17) 


In  eq.  (3.17),  £|^,  ^((?')  represent  the  components  of  the  free  space  dipole 
field  at  in  the  guiding  edge  fixed  coordinate  system. 
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3.5  An  equivalent  current  approach 

The  surface  radiation  integral  in  eq.  (3.8)  asymptotically  reduces  to  a  line 
integral  along  edges  (1)  and  (2)  of  the  trihedral  structure.  Thus,  the  high 
frequency  approximation  of  the  edge  and  vertex  diffracted  fields  can  be 
viewed  as  the  stationary  phase  and  end  point  contribution,  respectively, 
of  the  radiation  integral  of  equivalent  line  sources  “excited”  by  the  im¬ 
pinging  edge  wave.  As  far  as  the  edge  wave  current  effect  developed  in 
Section  3.4  is  concerned,  the  same  expression  for  the  field,  excluding  the 
factor  F\2kLc  sin’(^/2)],  can  be  obtained  via  the  asymptotic  approximation 
of  the  line  radiation  integral  of  the  physical  optics  components  of  Michaeli’s 
equiAralent  currents  along  edges  (1)  and  (2)  and  being  proportional  at  each 
point  to  the  tangential  components  of  the  incident  edge  wave.  This  ob¬ 
servation  suggests  the  addition  of  a  fringe  effect  to  the  edge  wave  current 
related  diffracted  field  which  is,  analogously,  described  by  the  radiation  in¬ 
tegral  of  Michaeli’s  fringe  currents.  Although  the  latter  have  been  derived 
for  plane  wave  incidence,  reasonable  results  are  obtained  if  one  generalises 
the  concept  for  arbitrary  (non-uniform)  ray  optical  wave  fronts  and,  hence, 
for  an  edge  wave. 

As  before,  our  attention  is  restricted  to  face  1  of  the  wedge.  The  equiv¬ 
alent  edge  currents  presented  in  [7]  for  grazing  edge  wave  incidence  at  the 
point  Q\  of  edge  (1)  may  be  expressed  in  the  form 


/.(«.) — 


■(Qi)U 


(3.18) 


«■((?.)  =  (3,19) 

t'l  is  the  rmit  vector  tangent  to  the  edge  at  Qi  and  ti,  mi  are  known  and 
slowly  varying  functions  of  the  observer’s  location  and  the  distance  of  Qx 
from  the  tip  kA  the  trihedron.  A  betor  of  1/2  has  been  also  introduced 
due  to  grazing  incidence. 

It  is  presumed  that  the  edge  and  vertex  diffracted  edge  wave  associated 
with  face  1  of  the  wedge  can  be  approximated  by  the  radiation  int^al  of 
the  equivalent  currents  Afi(Qi)  “flowing”  along  edge  (1),  which,  in 

the  FVesnel  or  the  Fraunhofer  region  of  the  edge,  can  be  explicitly  written 
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ais 


E^{s,)  «  jk''  Zo  sin*'  ai  C(i/)  At{p',  4>') 

•  /  [fi  X  X  + X  <im,(j;;<i)] 

Jo 

,u-x  exp{-jk{R  +  ri)} 

with  R  =  lla'c  —  ^1^*1  II  nnd  ^  =  (i*c  —  iiii)fR. 


(3.20) 


3.5.1  Edge  wave  vertex  diffracted  field 

The  phase  of  the  integrand  in  (3.20)  exhibits  a  stationary  phase  point  in  the 
neighborhood  of  an  end  point,  the  latter  being  coincident  with  the  branch 
singularity  of  the  integrand.  Therefore,  for  large  values  of  the  parameter  k 
the  asymptotic  evaluation  of  the  integral  reveals  a  vertex  contribution  as 
well  as  an  edge  diffracted  term,  if 

cos  oi  ( 1  —  cos  /9)  +  sin  ai  sin  cos  ^  >  0, 

i.e.,  if  Q*  lies  on  edge  (1)  itself  rather  than  on  its  extension.  From  the 
leading  term  in  the  asymptotic  expansion  of  the  integral  in  (3.20),  one 
obtains  the  following  approximation  of  the  vertex  diffracted  field  (associated 
with  face  1): 

El\s\)  %  exp{j(i/  +  l)’r/2}  r(i/)  C{u)  Zo  Al(p',  4>') 
exp{-jk{Bc  +  s^)}  G\{ac)  FlJ^ka\{a^\  s;)) 
iiraciKy  [cotai(l  —  cos/?) sin/? cos 

in  which  we  have  substituted 

Gj(sc)  =  ic  X  ie  X  <iti(Q')  +  ic  X  <imi((?‘').  (3.22) 

The  previous  discussion  also  suggests  the  decomposition 


J  i{(Q’tl‘{)  \. 

1  iriQ'iMD  1 

1  rnx{Q‘)  j 

1  ’"{(Q'i/'i)  / 

'  1  1 

*i(Q'),  ”»i(Q').  can  be  found  in  [7]  along 

with  the  definitions  (3.18),  (3.19). 
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The  parameters  involved  in  the  expressions  of  the  edge  equivalent 
currents  depend  on  the  choice  of  the  edge  fixed  coordinate  system  and,  in 
particular,  on  the  angle  between  the  tangent  to  the  edge  at  Qi  and  the 
unit  sector  &  tangent  to  face  1  at  Qi  which  may  be  chosen  arbitrarily. 
When  the  currents  are  truncated  to  reveal  an  end  point  effect,  one  should 
be  cautious  about  the  choice  of  the  parameters  ,  i.e.,  the  choice  of 
the  edge  fixed  coordii-ate  system,  so  that  the  terminated  equivalent  edge 
currents  represent  correctly  the  end  point  effect  of  the  truncated,  by  the 
edges  of  the  trihedron,  true  surface  induced  currents.  Within  this  context, 
a  correct  choice  of  the  parameter  for  the  physical  optics  component  of 
the  equivalent  currents  is 


= 


sin  tti  sin  /3i  cos  +  cos  oi  (cos  +  cos  Oi ) 


sin*  Oi 


(3.24) 


so  that  the  unit  vector  &  is  parallel  to  the  guiding  edge,  whereas  for  the 
fringe  currents  the  proper  choice  is 


sin  oi  sin  0i  cos  <f>i  —  cos  oi  (cos  +  cos  Oi ) 
sin*  oi 


(3.25) 


and  now  a  is  parallel  to  the  edge  diffracted  ray  from  Q‘  that  grazes  face  1. 

Not  surprisingly,  using  eq.  (3.24)  for  the  definition  of  the  parameter 
in  the  expressions  of  the  physical  optics  equivalent  edge  currents  one  finds 

(3.26) 


so  that  the  vertex  diffracted  field  assumes  the  approximation 

(3.27) 


in  which  Ei^{sc)  is  identical  with  the  term  associated  with  face  1  in  (3.17), 
while  the  fringe  current  related  wave  equals 

Ej*^(i;)  «  exp{i(i/+  l)w/2}r(t/)C'(i/)Zo34:(p',^')  sin^Qx 

cxp{-jk{sc  s;)}  F{2kLc sin*  f )  Ff  „[feo;(s;;  s^)] 

\lTBe{a'gY  (cOSQi  +  COS)3i)*' 

(3.28) 
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The  vector  function  Gi’^(se)  is  related  with  Michaeli’s  equivalent  currents 
at  Q‘  via  the  equation 


X  «c  X  x  ii. 


(3.29) 


The  above  expressions  simplify  considerably  for  the  case  of  the  plane 
angular  sector,  with  oi  =  02  =  a  ,ai  =  Oj  =  o.  Without  presenting  the 
details  of  the  derivation,  which  involves  only  elementary  manipulations,  for 
the  field  related  with  the  fringe  currents  one  finds 


G‘’^(sc;  a)  F(2kLc  sin’  f )  Ff  «c)] 

y/cot  a(l  —  cos  /3)  4-  sin  /3  cos  <f> 


where  the  function  G^‘^(3e)  is  given  by 

1 


G^’H^c)  = 


cos(2a)  +  cobS 


-/fsin^  +  + 


sin’ o  sin  ^1  cos(j/2) 


with 


cos  S  =  sin(2a)  sin  0  cos  ^  —  cos(2a)  cos  /?, 


(3.30) 


(3.31) 

(3.32) 


<t>)  =  —(cos  a  +  cos/9i)  (sin  a  cos cos  <f>  +  cos  a  ein/S) 
•[sin  a  cos  ^1(1  —  cosacosySi)  —  cos’asin^i] 

+  sin‘ a  sin  ^  sin  ,  (3.33) 


<!>)  =  sin  a  sin  <f>  (cos  a  -f  cos  0i ) 

•[sin  a  cos  ^1(1  —  cosacos^di)  —  cos’ a  sin /?i] 

+  sin^  a  sin  (sin  a  cos ^  cos  ^  -f  cos  a  sin  fi).  (3.34) 

For  the  particular  case  of  the  right  angle  sector  (a  =  7r/2),  it  readily 
follows  that 

C?'-'(i;;»r/2)  =4(cot|  -  csc^)  (3.35) 
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and  the  total  tip  diffracted  field  assumes  the  simple  representation 


-jZo  +  «e)} 

“  ■^17  ^ 

F(2kL^m‘  j)  f^i/i(tL,(l  -  ^1  -Bin’^cos*^)!  . 


^sin  a  cos  ( 


13.  (3.36) 


3.5.2  Edge  wave  edge  diffracted  field 

It  can  be  shown  that  the  edge  diffracted  field,  i.e.,  the  stationary  phase 
contribution  to  the  radiation  integral  of  Michaeli’s  equivalent  currents,  is 
the  edge  diffracted  edge  wave  predicted  by  UTD,  multiplicatively  corrected 
by  the  transition  function 

CO.’ 

where 

K-i()  =  (3-3T) 

with  the  star  denoting  complex  conjugate.  The  large  parameter  equals 

r  - 

The  transition  function  guarantees  the  finiteness  of  the  edge  dif¬ 

fracted  edge  wave  as  Qj  — »  and  the  unifonmiy  (continuity)  of  the  to¬ 
tal  field  as  the  “shadow  boundary'’  of  the  edge  diffracted  rays  (the  cone 
[cotai(l  -  COB0)  -I-  sin^Scos^]  =  0)  is  crossed.  It  appears  as  a  type  of 
caustic  correction  factor  in  the  sense  that  it  compensates  the  singularity  of 
the  edge  diffracted  field  at  the  extension  of  the  guiding  edge.  However,  in 
our  case  the  singularity  arises  from  the  behavior  of  the  incident  field  on  the 
edge  rather  than  the  focusing  of  the  diffracted  rays  into  a  caustic. 


(3.38) 


3.5.3  A  heuristic  correction  factor  for  the  direct  edge 
wave 

As  iJready  pointed  out,  the  rigorous  solution  of  the  radiation  of  a  dipole 
in  the  presence  of  an  infinite  wedge  predicts  a  singular  field  at  the  edge 
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of  the  wedge,  in  consistency  with  Meixner’s  edge  condition.  However,  for 
a  semi-infinite  structure  (trihedron),  the  singularity  of  the  edge  wave  at 
the  extension  of  the  guiding  edge  »  0)  in  free  space  is  not  physically 
acceptable.  To  overcome  this  discrepancy,  the  multiplicative  correction  of 
the  edge  wave  associated  with  a  semi-infinite  or  finite  edge  with  the  use 
of  a  suitable  transition  function  (in  general  a  dyadic)  has  been  suggested. 
Such  a  transition  function  can  be  empirically  derived  by  requiring  the  con¬ 
tinuity  of  the  total  field  at  the  shadow  boundaries  of  the  direct  wave  (edge 
wave),  namely  at  the  planes  This  continuity  was  guaranteed  by 

the  UTD  evaluation  of  the  edge  diffracted  field,  but  it  is  violated  in  the 
paraxial  region  after  the  introduction  of  the  function  which  assures 

the  uniformity  of  the  total  diffracted  field.  Obviously,  the  edge  wave  can 
be  midtiplied  by  a  similar  transition  function  so  that  the  total  field  re¬ 
tains  its  continuity  in  the  paraxial  region  as  well  as  outside  of  it,  where 
the  transition  function  approaches  unity.  In  addition,  such  a  multiplicative 
correction  would  yield  a  finite  total  field  along  the  extension  of  the  edge  in 
free  space. 

In  particular,  for  the  plane  angular  sector,  a  convenient  modification  of 
the  edge  wave,  in  the  extension  of  the  guiding  edge,  reads 


2iry/2v 


■00  CO.  A  .in(*/2)  +  i  com)]  2E(z2^. 

y  Sin  Po  So 


(3.39) 


with  So  denoting  the  distance  of  the  observation  point  from  the  dipole  and 
00  is  the  elevation  angle  of  the  observer  in  the  guiding  edge  fixed  coordinate 
system  centered  at  the  point  of  the  projection  of  the  dipole  onto  the  edge. 


3.6  Discussion  and  numerical  results 

The  edge  wave  vertex  diffruted  field  is  a  higher  order  term  with  respect  to 
the  large  parameter  k  in  the  asymptotic  solution  of  the  radiation  of  a  dipole 
in  the  vicinity  of  the  edge  of  a  trihedron.  Nonetheless,  it  contributes  signif¬ 
icantly  to  the  ^-directed  component  of  the  field,  especially  in  the  paraxial 
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region  of  the  guiding  edge  and  along  its  extension.  This  is  due  to  the  ac¬ 
cumulation  of  electric  current  flow  lines  in  the  vicinity  of  the  guiding  edge 
excited  by  the  ray  optical  edge  wave. 

The  solution,  based  on  Michaeli’s  equivalent  currents,  is  essentially  an 
asymptotic  PTD  (Physical  Theory  of  Diffraction)  approach,  in  that  a  fnnge 
current  effect  due  to  the  terminating  edges  is  added  to  the  edge  wave  cur¬ 
rents,  cast  finally  into  a  UTD  form.  Although,  the  derivation  of  the  fringe 
edge  currents  assume  an  infinite  edge  and  uniform  plane  wave  illumination, 
the  field  is  expected  to  retain  its  singular  behavior  in  the  vicinity  of  the 
vertex,  which,  moreover,  does  not  contradict  the  “tip  condition”  (i.e.,  the 
behavior  of  the  field  in  the  neighborhood  of  a  vertex)  as  investigated  by  sev¬ 
eral  authors.  It  should  be  emphasized  that  the  approximations  attempted 
in  Sections  3.4  and  3.5  by  no  means  present  a  complete  rigorous  represen¬ 
tation  of  the  tip  diffracted  field,  but  it  merely  includes  the  information  of 
the  truncation  of  known  components  of  the  currents  flowing  over  the  wedge 
surface. 

The  patterns  of  the  total  diffracted  field  predicted  from  the  edge  wave 
current  radiation  integral  and  the  equivalent  source  approximation,  for  the 
configuration  shown  in  Fig.  3.4(b),  are  compared  in  Fig.  3.4(a).  Clearly,  the 
tip  effect  of  the  edge  wave  currents  does  not  compensate  the  discontinuity 
of  the  edge  wave  edge  diffracted  field.  On  the  other  hand,  the  equivalent 
currents  result  yields  a  continuous  pattern  across  the  shadow  boundary  cone 
of  the  terminating  edge  diffracted  rays,  and  appears  as  a  more  complete 
representation  of  the  diffraction  effects. 

Comparisons  of  the  calculated  field  (which  includes  only  the  two  edge 
wave  tip  diffracted  terms  corresponding  to  the  two  adjacent  comers  added 
to  the  modified  direct  edge  wave)  via  Michaeli’s  equivalent  currents  (de¬ 
noted  as  MEC  on  the  graphs)  and  the  edge  wave  current  approach  (EWC) 
with  Method  of  Moments  results  (MM)  are  made  in  Figs.  3.5-3.7,  for  a 
4A  square  plate.  The  (fipole  has  been  placed  in  the  close  vicinity  of  one 
of  the  edges  (p'  =  O.OIA)  of  the  plate  and  sufficiently  far  from  its  comers 
(s'  =  2A).  The  direct  edge  wave  (EW)  (eq.  (3.3))  is  also  plotted  so  that 
the  effect  of  the  corner  diffracted  fields  is  better  illustrated.  Note  that  in 
Fig.  3.7  (where  ^  =  45*’)  the  edge  wave  edge  diffracted  fields  should  be  also 
added  to  the  comer  diffracted  fields.  The  agreement  is  good,  especially  in 
regions  where  contributions  from  other  diffraction  mechanisms  other  than 
the  comer  adjacent  to  the  guiding  edge  are  known  to  be  neglipble.  In  fact, 
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as  shown  in  Pigs.  3.5-3. 7,  the  agreement  between  the  Method  of  Moments 
data  and  the  calculated  field  progressively  improves  for  larger  values  of  the 
<j)  angle.  For  4>  =  3.7),  other  mechanisms  such  as  diffraction  from 

the  remote  corners  of  the  square  plate  as  well  as  double  and  triple  edge 
diffraction  may  contribute  significantly  to  the  pattern.  Our  approxima¬ 
tion  also  improves  when  the  paraxial  re^on  is  approached  (/9  — »  0°,  ISO**), 
where,  as  a  matter  of  fact,  the  total  field  is  stronger.  The  latter  justifies 
the  validity  of  the  asymptotic  analysis  which  resulted  in  the  multiplication 
of  two  transition  functions  as  well  as  the  choice  of  the  correction  factor  for 
the  direct  dipole  field.  Unfortunately,  the  ^-directed  field  does  not  exhibit 
an  analogous  agreement  (Fig.  3.8),  at  lower  azimuthal  cuts,  mainly  due 
to  the  fact  that  other  mechanisms  contribute  significantly  to  the  pattern 
of  that  polarization.  However,  the  small  angular  variation  of  the  pattern, 
which  is  a  typical  ^  component  pattern  for  ^  >  90*^,  indicates  that  the 
corner  associated  ^-directed  fields  are  sufficiently  weaker  compared  with 
the  ^  component  of  the  direct  edge  wave  and  can  be  neglected  in  practical 
calculations. 

The  resiilts  also  reveal  a  small  variation  of  the  total  field  with  respect 
to  the  azimuthal  coordinate  (angle  ^),  in  contrast  with  the  relative  large 
changes  of  the  calculated  pattern  in  the  elevation  plane.  The  <f>  dependence 
becomes  significant  only  at  lower  cuts  where  the  contribution  from  the 
opposite  edge  and  its  two  adjacent  comers  is  appreciable. 

The  second  example  examined  also  involves  the  radiation  of  a  small 
monopole  in  the  close  vicinity  of  one  of  the  edges  of  a  rectangular  plate 
(p'  =  O.OIA),  but  now  with  different  distances  from  the  two  adjacent  comers 
(s'j  =  2A,  s'j  =  3A).  The  calculated  field,  which  again  includes  two  tip 
diffracted  rays  in  addition  to  the  direct  edge  wave,  is  compared  with  Method 
of  Moments  results  as  well  as  measured  data,  as  shown  in  Figs.  3.9  and  3.10 
in  the  azimuthal  planes  ^  =  180”  and  ^  =  150”,  respectively.  The  accuracy 
of  the  measurement  deteriorated  at  smaller  azimuthal  angles,  where  the 
support  structure  influenced  significantly  the  measured  radiation  pattern. 
The  measured  data  in  Fig.  3.9  was  obtained  with  an  absorbing  material 
(RAM)  around  the  remote  corners  and  the  opposite  edge  of  the  rectangular 
plate,  so  that  the  pattern  is  more  closely  related  to  the  calculated  field, 
which  again  includes  two  adjacent  tip  diffracted  rays  in  addition  to  the 
direct  edge  wave.  Clearly,  the  agreement  is  better  in  this  case,  in  contrast 
with  the  results  of  Fig.  3.10  where  the  RAM  has  been  removed. 
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It  should  be  noted  that  neither  the  edge  wave  current  solution  nor  the 
equivalent  source  formulation  axe  expected  to  yield  accurate  results  for 
sm2dl  angular  sector  angles  (a  <<  ir/2).  In  this  case  a  strong  coupling 
between  the  two  edges  forming  the  sector  occurs,  which  is  not  encoimtered 
in  the  evaluation  of  the  radiation  integral  of  the  currents  flowing  along  the 
edge  (1).  Moreover,  when  a  — »  ir,  the  edge  wave  currents  associated  tip 
diffracted  field  vanishes  smd  the  total  solution  reduces  to  the  edge  wave  over 
an  infinite  wedge,  whereas  the  same  property  is  not  true  for  the  equivalent 
current  formulation  of  the  vertex  diffracted  field.  The  latter,  therefore,  fails 
in  cases  of  very  wide  angles,  which  require  a  more  careful  treatment. 

It  should  be  finally  noted  that  the  analysis  related  with  a  dipole  excited 
edge  wave  can  be  easily  extended  to  treat  the  diffraction  of  an  edge  wave  of 
the  form  (3.6).  However,  the  problem  of  a  plane  wave  at  grasing  incidence 
on  an  infinite  or  semi-infinite  wedge  is  of  rather  less  practical  importance 
than  that  of  edge-on  incidence  for  a  finite  structure. 


3.7  Summary 

Our  major  objective  was  to  describe  approximately  the  edge  wave  diffrac¬ 
tion  mechanisms  associated  with  the  interaction  of  an  edge  wave  and  the 
vertex  of  a  trihedron  as  a  first  step  towards  the  study  of  more  complex  ge¬ 
ometries.  To  approximate  the  vertex  diffracted  wave  a  UTD  solution  was 
developed  based  on  a  PTD  like  approach  to  defining  the  currents.  The  addi¬ 
tion  of  a  fringe  current  component  flowing  in  the  vicinity  of  the  terminating 
edges  yields  a  continuous  radiation  pattern  across  the  shadow  boundaries 
of  the  edge  diffracted  waves.  Although,  the  approach  is  neither  rigorous 
nor  complete  from  a  PTD  point  of  view  it  yields  comparable  results  with 
the  Method  of  Moments  as  well  measured  data  and  can  be  used  as  a  first 
approximation  to  the  edge  wave  diffraction  problem. 

The  extension  of  the  results  of  this  Chapter  for  the  case  of  tip  excited 
edge  waves  and  the  study  of  the  edge  wave  diffraction  mechanisms  associ¬ 
ated  with  an  edged  body  illuminated  by  a  plane  wave,  which  are  of  partic¬ 
ular  importance  when  edge-on  incidence  is  approached,  is  investigated  in 
Chapter  4. 
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90.0 


(a) 


Figure  3.4:  (a)  Angular  section  Geometry,  (b)  Far  region  edge  and  vertex 
diffracted  /^-directed  field  at  the  cone  0  =  85°:  dotted  line  -  edge  wave 
current  solution,  solid  line  -  edge  wave  and  fringe  current  solution. 
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BETA  ANGLE  (DEG.) 


Figure  3.5:  Far  field  /9-directed  wave  at  the  asimuthal  cut  ^  =  180®  (nor 
malized  to  2.3dB). 


BETA  angle  (DEG.) 


Figure  3.6:  Far  field  /9-directed  wave  at  the  azimuthal  cut  ^  =  135®  ( 
malized  to  1.9dB). 


BETA  angle  (DEG.) 


Figure  3.8:  Far  field  ^directed  ware  at  the  azimuthal  cut  ^  =  120”. 


BETA  ANGLE  (DEG.) 


Figure  3.9:  Far  field  /^-directed  wave  at  the  asimuthal  cut  ^  =  180°  ( 
malised  to  5dB). 


Chapter  4 

Edge  Wave  Mechanisms  for 
Polyhedral  Structures 

4.1  Introduction 

Experimental  evidence  (e.g.,  impulse  response  of  a  square  plate)  indicates 
that  edge  wave  diffraction  mechanisms  become  significant  when  edge  on 
incidence  is  approached  and,  in  many  circumstances,  dominate  the  lower 
level  regions  of  the  monostatic  and  bistatic  cross  section  of  edged  struc¬ 
tures.  Therefore,  the  incorporation  of  this  type  of  mechanisms  in  cross 
section  computations  serve  as  an  additive  correction  to  the  comer  diffrac¬ 
tion  solution  developed  by  Brinkley  and  Marhefka  [2]. 

The  edge  wave  diffraction  phenomenon  on  flat  plates  has  been  recently 
studied  by  Sikta  [11].  His  solution,  although  it  provided  reasonable  agree¬ 
ment  in  several  cases,  is  based,  as  already  pointed  out,  on  an  empirical  cor¬ 
ner  diffraction  expression,  while  his  study  of  multiple  interactions  reqtiired 
the  heuristic  introduction  of  a  ‘^reflection  coefBcient”.  A  brief  outline  of  his 
formulation  of  the  problem  is  presented  in  Section  4.2,  so  that  comparisons 
against  the  present  development  can  be  made. 

In  Section  4.3  the  tip  excited  edge  wave  and  the  associated  interactions 
between  two  adjacent  vertices  of  a  polyhedral  structure  are  examined,  em¬ 
ploying  the  more  complete  depiction  of  the  corner  effect  developed  in  Chap¬ 
ter  3,  along  with  reciprocity,  and  based  on  plausible  heuristic  arguments. 
Specifically,  in  high  frequencies  the  vertex  formed  by  the  intersection  of 
three  plates  can  be  considered  as  the  tip  of  an  ideal  trihedron  (semi-infinite 
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wedge).  Hence,  when  a  ray  optical  field  impinges  onto  the  tip,  an  edge 
wave  is  excited  that  can  be  determined  via  reciprocity  &om  the  expressions 
of  the  edge  wave  vertex  diffracted  field.  The  edge  wave  itself  is  ray  optical 
and  hitting  an  adjacent  corner  it  is  rediffracted  accordingly. 


Figure  4.1:  Comer  excited  edge  wxvt. 


4.2  Sikta’s  formulation 


The  edge  wave  excited  by  the  tip  and  propagating  along  the  edge  of 
the  plate  is  as  shown  in  Fig  4.1.  According  to  Sikta,  it  is  determined  to  be 

I  1  I  Bj,.  ]  ,QZexp(-jfcn  (4 « 

I  fi.  on*'.; j  V '  ' 

in  which  the  edge  wave  launching  coefficients  are  defined  by 


Dr{B'c;iM  = 

h 


exp(-jw/4)  1  -1-  cob/9' 


2\/2^  nn{/3'/2) 

j  sin(^72)sin(^./2)  \ 
!/  \  co8(^72)  co«(^./2)  J 


(4.2) 
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In  the  above  expressions 


(4.3) 


s' ,  j9',  <!>'  are  the  coordinates  of  a  dipole  source  in  the  edge  fixed  coordinate 
centered  at  ,  while  I  is  the  length  along  the  edge  from  the  tip  and  tf>c  is 
the  azimuth  smgle  of  a  point  close  to  the  edge.  Note  that  Sikta’s  expression, 
although  it  appears  to  have  the  correct  <f>c  and  I  dependence  characteristic 
of  an  edge  wave,  it  does  not  exhibit  the  proper  edge  singularity. 


Figure  4.2:  Geometry  of  the  corner-to-comer  edge  wave  interaction. 

Two  adjacent  comers  of  a  flat  plate,  along  with  the  involved  geometric 
quantities,  su'e  illustrated  in  Fig.  4.2.  The  source  point  is  again  located 
at  with  respect  to  the  guiding  edge  fixed  spherical  co-.  .dinate 

system  centered  at  the  tip  while  the  location  i2(sei/9i  4")  of  receiver 
is  defined  in  terms  of  an  edge  fixed  coordinate  system  with  its  origin  at 
the  point  The  field  .E'*®*,  produced  by  edge  wave  diffraction  at  Q®*, 
the  edge  wave  being  excited  by  the  tip  due  to  source  illumination  and 
guided  along  the  edge,  is  of  interest  here.  It  is  noted  that  the  total  double 
corner  interaction  should  include  also  a  ray  along  the  path  SQ'^Q^'R. 

The  far  field  analysis  of  the  edge  wave  mechanisms,  always  in  the  con¬ 
text  of  Sikta’s  analysis,  involves  the  multiplication  of  the  incident  field 
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Ep,,  E^i  with  the  diffraction  coefficients  in  eq.  (4.?)  and  a  reflection  coeffi¬ 
cient,  which,  for  a  right  amgled  corner,  should  read 


exp(-jir/4) 


(4.4) 


Specifically,  the  edge  wave  diffraction  by  two  adjacent  comers  in  a  rectan¬ 
gular  plate  becomes 


EfiW^)  \  +  <c)] 


\  E^{Q‘^)  I 


(4.5) 


Explicitly, 


Ep^^  \  _  -  exp(jTr/4)  exp(-iA;(f  +  Je) 

El^‘^  f  (2^ib)>/» 

cos»()9/2)F[2ibf  An\fil2)\  co8>(/g72)F[2fef  8in»(/972)] 
sin(^/2)  8in(/?V2) 

/  sin(^/2)sin(<^72)  \/28in(^/2)cos(^72)  1  f  E^>(Q^')  1 
(  v^co8(^/2)8in(^72)  2co8(0/2)co8(072)  /  \  E^,(Q‘')  J' 

(4.6) 


4.3  Edge  wave  interactions  at  a  polyhedral 
structure 

Assuming  that  the  point  Q*’*  is  the  tip  of  a  semi-infinite  wec^e,  which 
is  justifiable  as  long  as  diffraction  is  a  local  phenomenon  (i.e.,  when  the 
length  of  the  edge  is  sufficiently  large  in  terms  of  the  wavelength),  the 
field  excited  by  this  comer  and  propagating  along  edge  (1)  can  be  obtained 
via  reciprocity  from  eq.  (3.27)  along  with  eqs.  (3.17)  and  (3.28).  Hence, 
considering  a  point  electric  source  with  moment  ^  located  at  the  close 
neighborhood  of  ,  then  it  has  been  shown  that  the  edge  wave  this  dipole 
excites  is  diffracted  by  the  tip  and  is  evaluated  at  5  according  to 

=  ZoexpL;(t/  -I-  l)ir/2]  r(*/)  C{t^) 
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p*  8in(^«/2)  -f  co8(^,/2) 

Pi-^ 

^  (4.7) 


where  (pe^^c)  O’^e  the  cylindrical  coordinates  of  the  point  source  at  the 


Figure  4.3:  Edge  wave  diffraction  for  a  polyhedral  structure, 
guiding  edge  associated  coordinate  system,  while 

G(b'c]  /;  Oei ,  n.i ;  0/1 ,  n/i )  =  F(2fcl'  sin*(/9/2)] 

sin^  Q,iG;,(s'c)F!^{fe£;[l  -  cos(aei  -  j9*i)]} 

(cos  0,1  —  COS^'j)*' 

I  »in‘^a/ig/i(Pe)F!^{fe£;[l  +  cos(a/i  + 

(cOSQ/i  +  COS/?^j)*'  *  ' 

and 

5:., /,(<".)  =  -  coi(;872)/f- + (4.9) 

with  Gli  ji  defined  by  eq.  (3.29),  adapted  to  the  coordinate  systems  related 
to  edges  (el)  and  (fl).  It  is  reminded  that  i/  =  1/n,  with  n  denoting  the 
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guiding  wedge  number.  Further,  n^\j\  and,  later,  n/3,.3  are  used  to  denote 
the  wedge  number  associated  with  the  a4jacent  edges  (el),  (fl)  and  (e2), 
({2),  respectively.  The  parameter  X'  is  defined  by  eq.  (4.3).  It  should  be 
finally  noted  that  the  angles  are  the  elevation  angle  of  S  in  edges  (el) 
and  (fl)  fixed  coordinate  systems  as  shown  in  Fig.  4.3. 

Considering  a  dipole  source  with  moment  ^  located  at 

S  and  applying  the  reciprocity  theorem  the  corner  excited  edge  wave  at  the 
location  of  the  point  source  equals 


E-  =  + i)x/2|  r(..)c(.) 

Pi 

exp[-jfc(s*  +  /)]  ^  p  ^  \\  r? 

-  exp(yi/ir/2)r(i/)C7(i/)  p*  sin(i/^e)  +  cos(i/^,) 

k  pl- 

.J‘(Q«).  (4,10) 


It  is  reminded  that  the  paraxial  edge  wave  excited  by  a  point  source  radi¬ 
ating  in  the  vicinity  of  the  edge  of  a  half  plane  and  evaluated  at  the  point 
(Pe,^e,0  ibl  >>  1  can  be  approximated  by 


p,  8in(l/^e)  +  C08(t/^e)  exp(-jfef) 


pI~'' 


(4.11) 


where  is  a  source  related  factor  in  view  of  eq.  (3.5).  Upon  comparison 
of  eqs.  (4.10)  and  (4.11)  one  readily  obtains 


i4**  =  i4*’(s'e;l;a.i,n,i;a/i,n/i) 

jYo exp(ii/Tr/2)r(»/) ^ 

—  ^  (*  C>  *1  ®«1>***1  >  ® '1 » **/i)  * "  ( )» 


(4.12) 


*  thus,  defining  a  fictitious  equivalent  point  source  located  at  the  tip 
and  exciting  the  edge  wave  given  by  eq.  (4.10).  One  then  considers  the 
diffraction  of  the  edge  wave  produced  by  the  equivalent  point  source  at  the 
corner  Cj®*.  Referring  to  eq.  (3.21)  one  writes 

F'*'*  =  Zo  exp[i(i/  l)w/2)  r(v)  Ciu)  A*’ 
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47ra'^^~'  (4*13) 

Upon  incorporation  of  eq.  (4.12)  into  the  above  relation  one  derives  the  the 
following  expression  for  the  edge  wave  diffracted  field: 

=  g“(y., I, oi. a,)  ■  ^(ly  )  ,  (4.14) 

where  the  dyadic  diffraction  coefficient  is  given  by 

D  (.„.,;J.a.,a,)= - - 

(?*=(«;;/;  a,  j ,  n,, ;  a/a,  n/i)^'(Pe;  h  oiei ,  n.i ;  a/i ,  ny  i ).  (4.15) 

In  the  case  of  a  flat  plate  we  have 

^  (^e,  a*j,  Tl*j;  Qyjjflyj)  =  20o(SeJ QE]) 

=  2  [- cotW2)/j  +  (?'•'(*-)] , 

<?*(sV,^;aei,n,i;a/i,nyi)  =  2G5(s*'e;/,ai) 

=  2(-cot(^72)/9'  +  <3'’'(.‘.)], 

where  is  defined  by  eq.  (3.31)-(3.34),  and 

C(.)  =  C(V2)  =  ^^. 

Hence,  eq.  (4.15)  becomes 

5"(P., a„a.)  =  {4.16) 

The  above  can  be  further  simplified  in  the  case  of  two  adjacent  right  angled 
vertices.  In  particular,  the  comer-to-comer  edge  wave  interaction  becomes 

ncc  exp(jw/4)  ,  exp(-;fc(/  +  s^) 

-l2^dfc)vrW  ) 

F[2klnv?{Pl2)]  Ft^i^\2kl{\  -  y/l  -  sin*jgcoB»^] 
sin(/?/2)  v^sin  P  cos  ^ 

F[2A!f8in^(^72)]  Ff^y,(2fel(l  -  -  sin*  cos* 

nn{0'/2)  y/nn/3*coB^'  ' 

(4.17) 
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Note  that  eq.  (4.17)  predicts  a  dominant  $  directed  field  which  is  physi¬ 
cally  justifiable,  since  the  edge  wave  currents  are  accumulated  in  the  vicin¬ 
ity  of  the  guiding  edge.  The  same  is  not  true  with  the  corrected  Sikta’s 
expression  (4.6).  The  two  equations  indicate  also  a  different  azimuth  depen¬ 
dence  and  a  different  elevation  angle  variation  approximately  by  a  factor 
of  v^cBC^'cBCjS.  As  a  matter  of  fact,  eq.  (4.17)  is  expected  to  predict  a 
stronger  field  when  edge  on  incidence  is  approached.  FHirther,  the  transi¬ 
tional  behavior  of  the  two  expressions  is  different. 

It  should  be  reminded  that  the  comer  diffracted  edge  wave  is  discon¬ 
tinuous  at  the  shadow  boundary  of  the  edge  diffracted  edge  wave  in  such 
a  manner  that  the  total  field  remains  continuous.  In  the  same  context,  in 
addition  to  the  tip-to-tip  edge  wave  diffraction  demonstrated  above,  tip- 
to-edge  and,  vice  versa,  edge-to-tip  diffracted  ray  contributions  may  occur 
at  certain  aspects  of  incidence  and  observation.  The  discontinuity  of  these 
rays  when  the  tip  and  the  edge  diffraction  point  coalesce  should  be  normally 
compensated  by  an  analogous  discontinuous  behavior  of  double  corner  dif¬ 
fracted  edge  wave.  The  analysis  pertinent  to  this  type  of  rays  is  similar  and 
it  is  not  repeated  here.  It  simply  involves  the  multiplication  of  the  edge 
wave  in  eq.  (4.11),  excited  by  the  equivalent  point  source  by  the  proper  edge 
wave  edge  diffraction  coefficient  in  view  of  the  development  in  Section  3.5. 


4.4  Discussion  and  numerical  results 

The  diffraction  by  two  adjacent  comers  of  a  polyhedral  stmctme  has  been 
studied  as  a  second  order  additive  correction  to  the  corner  diffracted  fields 
in  a  high  frequency  presentation  of  the  scattering  by  conducting  edged  bod¬ 
ies.  A  double  coraer-to-corner  edge  wave  dyadic  coefficient  was  deduced, 
which  can  be  used  as  a  first  engineering  approximation  in  higher  order 
UTD  computations.  It  involves  the  multiplication  of  two  transition  factors 
expressed  in  terms  of  the  parabolic  cylinder  function  of  order  —v  and  the 
ordinary  FVesnel  integral,  that  appear  in  the  first  order  theory,  in  an  analo¬ 
gous  manner  that  the  douUe  edge  diffraction  is  described  via  ordinary  UTD 
calculations.  Outside  its  transition  domains  (i.e.,  when  the  source  or  the 
observer  are  lined  up  with  the  edge),  the  edge  wave  diffractkm  coefficient 
is  proportional  to  while  it  decays  with  the  inverse  distance  from  the 
tip  along  the  edge.  Therefore,  it  represents  one  of  the  next  higher  order 
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mechanisms  in  a  high  frequency  depiction  of  the  fields,  in  addition  to  the 
corner  diffracted  rays. 

The  effect  of  the  introduction  of  edge  waves  as  a  higher  order  mecha¬ 
nism  has  been  evaluated  via  RCS  computations  of  flat  plates.  The  comer 
diffraction  is  represented  as  the  truncation  effect  of  the  Physical  Optics 
and  non-uniform  component  of  the  currents  on  the  plate.  Two  flat  plate 
examples,  where  the  edge  wave  contribution  have  been  expected  to  be  of 
particular  importance  have  been  investigated:  (i)  the  off  principal  plane 
RCS  of  a  2A  square  plate  and  (ii)  Sikta’s  triangle.  Specifically,  Figs.  4.6- 
4.9  demonstrate  the  effect  of  the  superposition  of  edge  wave  mechanisms 
in  the  echo  of  a  2A  square  plate.  For  comparison,  the  corner  diffraction 
solution  [2]  is  also  shown. 

Sikta’s  triangle  (Fig.  4.4  and  4.5)  is  another  case  where  edge  waves 
become  important,  especially,  at  nose-on  incidence.  Indeed,  incorporation 
of  the  edge  wave  mechanisms  appears  to  improve  the  calculated  RCS  at 
regions  where  they  are  expected  to  be  significant.  It  should  be  emphasized 
that  corner-to-edge  (and  vice  versa)  edge  wave  mechanisms  which  may 
be  important  in  certain  aspects  of  observation  have  not  incorporated  in 
the  calculations.  Therefore,  a  discontinuity  may  appear  at  observation 
directions  corresponding  to  the  shadow  boundaries  of  the  corner-to-edge 
diffracted  rays.  Further  research  should  investigate  these  mechanisms  and 
their  transitional  behavior  as  well  as  their  relation  with  the  rays  arising  from 
the  asymptotic  integration  of  the  non-uniform  components  of  the  currents 
studied  in  Chapter  5. 

In  general,  due  to  the  non-availability  of  a  rigorously  derived  comer 
diffraction  coefficient,  one  is  difficult  to  assess  the  contribution  of  higher  or¬ 
der  mechanisms  from  RCS  pattern  studies.  Therefore,  time-domain  mech¬ 
anism  extraction  has  been  investigated  [37]  in  the  identification  and  quan¬ 
titative  analysis  of  the  doubly  diffracted  ray  contributions  that  constitute 
the  echo  of  a  large  edged  body.  It  has  thus  been  observed  that  the  edge 
waves,  at  least  outside  the  transition  regions  where  more  than  one  mech¬ 
anisms  coalesce,  are  adequately  predicted  with  the  present  solution.  In 
transitions  regions,  e.g.  at  edge  on  incidence  or  at  the  shadow  boundaries 
of  comer-to-edge  and  edge-to-comer  diffracted  rays  a  more  careful  investi¬ 
gation  is  required.  Fig.  4.10  depicts  the  measured  and  calculated  impulse 
response  of  a  diamond  shape  plate.  FVirther  investigation  of  the  edge  wave 
mechanisms  has  been  focused  in  the  evaluation  of  the  several  mechanisms 
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appearing  in  the  impulse  response  of  a  24"  square  plate  for  several  aspects 
of  illumination,  which  are  illustrated  in  Figs.  4.11-4.25.  The  comer  dif¬ 
fracted  rays  are  denoted  by  C1-C4,  while  possible  edge  wave  mechanisms 
are  indicated  by  E1-E4  and  double  diffractions  with  D1-D4.  The  support 
structure  returns  are  SI  and  S2. 
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Figure  4.6:  ff##*  RCS  of  a  2A  square  {date  at  •  =  30®. 


Figure  4.7:  RCS  of  a  2A  square  plate  at  ^  =  30®. 
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Figure  4.8:  RCS  of  a  2A  square  plate  at  =  45®. 


Figure  4.9:  o’M'  RCS  of  a  2A  square  plate  at  0  =  60®. 
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Figure  4.10:  Impulse  response  of  a  diamond  shape  plate  at  nose-on:  (a) 
measured,  (b)  calculated. 
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Figure  4.11:  Impulse  response  of  a  24”  square  plate  at  ^  =  45°  and  0  =  15°. 
Vertical  polarization,  (a)  Measured,  (b)  Calculated. 
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Figure  4.12:  Impulse  response  of  a  24”  square  plate  at  ^  =  45®  and  9  =  30®. 
Vertical  polarization,  (a)  Measured,  (b)  Calctilated. 
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Figure  4.16:  Impulse  response  of  a  24”  square  plate  at  ^  =  45°  and  9 
Horieontal  polarization,  (a)  Measured,  (b)  Calculated. 
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Figure  4.18:  Impulse  response  of  a  24”  square  plate  at  ^  =  45®  and  B  =  45°. 
Horizontal  polarization,  (a)  Measured,  (b)  CalcuSted. 
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Figure  4.20:  Impulse  response  of  a  24”  square  plate  at  ^  =  45®  and  =  75®. 
Horizontal  polarization,  (a)  Measured,  (b)  Calculated. 
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Figure  4.21:  Impulse  response  of  a  24”  square  plate  at  0  =  45®  and  9  =  15®. 
CrosB  pol (a)  Measured,  (b)  Calculated. 
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Figure  4.22:  Impulse  response  of  a  24”  square  plate  at  ^  =  45**  and  9  =  30*’. 
Cross  polarization,  (a)  Measured,  (b)  Calculated. 
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Figure  4.25:  Impulse  response  of  a  24”  square  plate  at  ^  =  45®  and  9  =  75®. 
Cross  polarization,  (a)  Measured,  (b)  Oalciilated. 
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Chapter  5 

A  Uniform  Ray  Approximation 
for  the  Interaction  Between 
Edges  of  Three-Dimensional 
Polyhedral  Structures 


5.1  Introduction 

A  uniform  ray  representation  of  the  interaction  between  edges  of  the  three- 
dimensional  far  field  region  scattering  by  a  flat  plate  structure  is  investi¬ 
gated.  This  is  accomplished  by  (i)  postulating  an  approximation  of  the 
induced  surface  ctirrent  on  ecM:h  face  composing  the  object,  and  (ii)  inte¬ 
grating  the  approximate  current,  either  in  closed  form,  or  asymptotically 
in  terms  of  the  well  tabulated  edge  transition  function.  Specifically,  the 
surface  current  on  each  plate  composing  the  structure  is  approximated,  in 
addition  to  the  usual  Physical  Optics  surface  current,  by  a  non-uniform 
current  excited  by  each  edge  of  the  plate,  which  is  derived  from  the  canon¬ 
ical  solution  (integral  representation)  to  the  infinite  wedge  problem,  and 
truncated  at  edges  of  the  plate.  The  superposition  of  a  secondary  non- 
uniform  current,  which  can  be  asymptotically  represented  by  a  Michaeli’s 
fringe  equivalent  source,  excited  by  doubly  diffracted  rays,  is  also  studied 
as  an  additive  correction  to  the  first  order  solution. 
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5.2  A  brief  review  of  Michaeli’s  equivalent 
currents 

A  plane  wave  illuminates  the  infinite  conducting  wedge  of  Fig.  5.1.  For 
convenience,  we  restrict  our  attention  to  face  1  of  the  wedge.  The  total 
current  is  the  superposition  of  the  well  known  Physical  Optics  current 

=  2n  X  J?‘(a:,y,z)|y=o  (5.1) 

and  a  non-uniform  component,  which  is  in  essence  the  current  excited  by 


n7T 


(7 

FACE  2 

Figure  5.1:  Infinite  wedge  illuminated  by  a  plane  wave. 


the  diffracted  fields;  it  is  depicted  in  Fig.  5.2  and  explicitly  described  by 
the  equations 


jHlf^cxpj-jkzcosfi') 

4irn 

.±1  ^di±lzll£ 


cxp(- jkx  sin/?' co8^)d^. 


j  exp(—jkz  cos 
4irn 
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(5.3) 


'2  /  [(-l)*’l^-E^o«C8c)9'8in^ +  if‘,cot/9'co8^] 

p.g=l 

•cot  ^2ti^  ^  }  exp(— jJb*8in/3'co8^)<i^. 

The  various  geometric  parameters  involved  in  the  above  expressions  are 


Figure  5.2:  Non-uniform  component  of  the  current  at  the  wedge. 


shown  in  Fig.  5.1.  The  contours  F,,  q  =  1,2,  are  the  steepest  descent 
SDP{±ir)  paths  through  the  scKidle  points  at  (  =  ±Tr  (Fig.  5.3).  Eqs.  (5.2), 
(5.3)  in  the  case  of  the  half  plane,  reduce  to  the  expressions 

>"“(x,z)  =  -^^^l4^^exp{j*(i8in^'co8(/>' -  zcos/3')] 

^  TT  81D  fj 

■{(yES^sin^'-^*,  cos  P'  cot  4>) 

I -  6'  4>' 

•F_  ( y  2fcz  sin  /?'|  cos  ^  |)  sign  (cos  ) 


"Tiro  f) 

•exp(-2jfci  sin cos*  y)|  > 


(5.4) 
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exp[7A;(®  sin^fl'cos^'  —  zcosp')] 

^/1r 

‘F-{<j2kx  sin  '0'\  cos  ^ |)  sign(coB  ^ ),  (5.5) 

iS  Ji 

that  involve  the  ordinary  Fresnel  integral 

F-(x)  =  j  exp(-jf<*)dt.  (5.6) 


Figure  5.3:  Steepest  descent  paths  in  the  complex  (  plane  (from  [Pathak 
and  Kouyoumjian,  1971]). 

Accordingly,  one  relates  a  non-uniform  field  with  the  above  currents. 
Provided  that  the  smallest  distance  of  the  field  point  from  the  surface  of  the 
wedge  is  sufRciently  large,  the  field  associated  with  j""  can  be  approximated 
by 

F""  »jkZ  I  RxRx  (5.7) 

Jsi  4itR 

with  Si  denoting  the  area  of  face  1  of  the  wedge.  It  has  been  shown  by 
Michaeli  [18,7]  that  the  above  surface  integral  can  be  reduced,  via  a  linear 
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phase  approximation,  to  a  line  integral  of  line  magnetic  and  electric  equiv¬ 
alent  sources  along  the  edge  of  the  wedge.  If  one  introduces  the  oblique 
system  of  coordinates  related  with  the  diffracted  rays  that  graze  the  wedge 
surface,  the  electric  and  magnetic  fringe  edge  currents  evaluated  at  the 
point  Q  of  the  edge  equal 


and 


(6.9) 


respectively,  where 


V{n  - 


-I-  COB^ 

(5.10) 


1  /  cot  ^  cos  ^ -I- cot  cos 

_  i - ..  .  - U{n  -  ) 


sin  p'  p  +  cos  ^ 

p cot /?'  —  cot P cos —  {l  —  4>' 


-I-  cot 


(^^)li-S3.  »■»> 


U{^  -  <!>•) 

p  -I-  cos^' 


sin  ^  J  I 
sin  0  sin  0'  \  | 

'  h(^^) 


2n  sin  7 


-f  cot 


The  expressions  for  the  edge  currents  pertinent  to  the  other  face  of  the 
wedge  are  obtained  by  means  of  a  simple  transformation.  The  parameter 
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(I  equals 


(6.13) 


/i  = 


cos  —  cos*  P' 
sin*/9'  ’ 


with 


cosV'  =  s'*  •  O’  =  sin sin )9 cos ^  +  cos )d' cos/?, 

and  O’  is  the  unit  vector  in  the  direction  of  the  grazing  diffracted  rays, 
whereas  the  angle  7  is  defined  by 


/X  =  cos  7, 


(5.15) 


along  with  the  (arbitrary)  choice  of  the  branch,  e.g.,  of  that  depicted  in 
Fig.  5.4.  This  particular  branch  selection  yields 


with 


7  =  -iln(^-f- 

-1) 

(5.16) 

a  1  _  / 

\  jVI  - 

A  V 

(5.17) 

It  becomes  clear  from  eqs.  (5.8),  (5.9)  and  eqs.  (5.10)-(5.12)  that  the 
only  circumstance  under  which  the  non-uniform  edge  equivalent  currents 
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become  singular  is  when  ^  =  0  (or  /*  =  1),  which  imply  that  v  =  i. 
This  singularity  is  integrable,  unless  simultaneously  =  w,  which  occurs 
at  glancing  incidence  and  forward  scattering  observation  (Ufimtsev’s  sin¬ 
gularity  [14]).  It  should  be  noted,  however,  that  Ufimtsev’s  singularity  is 
eliminated  when  one  integrates  the  non-uniform  currents  along  a  finite  path 
in  the  a  direction. 

Michaeli’s  total  currents  on  the  other  hand  exhibit  discontinuities  as  a 
result  of  the  non-uniform  asymptotic  approximation  of  the  physical  optics 
surface  currents.  In  general,  for  an  infinite  wedge,  it  is  impossible  to  obtain 
uniform  total  edge  currents  if  linear  phase  approximation  is  introduced 
for  the  reduction  of  the  surface  radiation  integral  to  a  line  integral  along 
the  edges.  Therefore,  for  the  infinite  wedge  case  one  has  to  assume  finite 
distances  of  the  observation  point  from  the  edge  and  introduce  a  quadrature 
phase  approximation,  which  will  yield  the  expressions 


2m^  -  <!>•) 

k  sin’  +  cos 


-  ff‘(Q)Mn^'(cot/3co6^  +  cot^'co,^')j 


ks  sin^  -I-  cos 
2  sin’  ^ 


(5.18) 


M^{Q)  =  - 


2jZHi{Q)U{n  - 
k  sin  P'  sin  -I-  cos  ^') 


F 


ks  sin^  P'{fM  -I-  cos  ^')* 
2  sin’  ^ 


(5.19) 


where  F(-)  is  the  familiar  edge  transition  function.  Besides,  in  case  of 
finite  bodies,  one  may  integrate  directly  and  effidently  the  Physical  Optics 
surface  currents  [38]. 


5.3  Higher  order  terms  using  non-uniform 
surface  wedge  currents 

The  total  surface  current  on  each  plate  of  the  structure  is  approximated,  as 
before,  by  the  Physical  Optics  current  and  a  non-uniform  current  excited  by 
each  edge  of  a  given  plate.  The  approximation  is,  in  some  respects,  similar 


78 


with  Michaeli’s  secondary  equivalent  current  development  [39].  However, 
the  latter  is  restricted  in  the  definition  of  a  secondary  edge  current,  which 
is  then  integrated  numerically  and,  in  addition,  no  explicit  expression  for 
the  corner  effect  is  given.  Here,  an  attempt  is  made  to  present  the  total 
solution  in  terms  of  ray  contributions,  corresponding  to  tip,  tip-to-edge  and 
edge  -to-tip  diffracted  rays.  These  rays  result  from  the  asymptotic  surface 
current  integration  and,  in  general,  do  not  satisfy  Fermat’s  principle  and, 
further,  radiate  as  if  the  body  was  not  present  and,  therefore,  they  are  not 
shadowed.  Nevertheless,  the  final  solution  is  expected  to  be  as  efficient  as 
GTD  in  high  frequency  calculations. 


Figure  5.5:  Non-uniform  current  excited  by  an  edge  and  associated  ray 
contributions  at  a  polygonal  structure  illiuninated  by  a  plane  wave. 

Referring  to  Fig.  5.5,  the  exact  non-uniform  surface  currents  excited  by 
edge  (1)  are  truncated  by  edge  (2),  provided  that  at  least  part  of  edge  (2) 
is  illuminated  by  the  diffracted  rays  emanating  from  edge  (1)  that  graze 
the  surface  of  the  plate.  Our  analysis  is  restricted  to  the  study  of  the 
non-uniform  current  emanating  from  edge  (1)  of  the  illustrated  face  of  the 
object;  similar  arguments  apply  for  any  other  edge  of  this  plate  and  for 
any  other  plate  of  the  body.  The  possible  existence  of  shadow  boundaries 
of  the  non-uniform  surface  current  associated  with  each  edge  should  be 
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emphasized.  The  discontinuities  of  the  surface  non-uniform  current  should 
be  normally  (in  a  PTD  sense)  compensated  by  an  analogous  discontinuous 
behavior  of  a  tip  excited  current,  which,  however,  remains  unknown  and  it 
is  not  included. 

One  establishes  the  oblique  coordinate  system  with  d  p<nnting 

to  the  direction  of  the  grazing  diffracted  rays  firom  edge  (1)  and  f  being 
parallel  to  the  edge  vector.  It  is  related  with  the  edge  (1)  fixed  cartesian 
system  of  coordinates  via  the  linear  transformation 


*1  =  <rsin/3|,  Zi  =  t  +  acoa0[.  (5.20) 

Note  that  subscript  1  in  the  angles  should  imply  reference  to  edge  (1)  fixed 
coordinate  system.  The  field  pertinent  again  to  the  non-uniform  current 
emsmating  from  edge  (1)  over  the  shown  face  of  the  object,  and  consistent 
with  our  postulate  for  the  total  current,  is  written  as 

ss  iibZsin/3{--^— X  s**  X 

r  j^(<T\r')txp{jkS^  •  P)da^dT\  (5.21) 
Jri  Jo 

since,  in  the  far  field, 

Rss  3  —  3^  '  {a' a  -I-  r'f)  =  s  —  a  con  if}  —  r  confix,  (5.22) 


with  3  being  the  distance  from  a  fixed  reference  point  on  edge  (1).  Hence, 
one  obtains 


jkZ  aTlfi[cXipi-jk3) ,  _  .  _  .  a  r\ 

= - — — - : - ^(cos^i  cos^i  /,  -  nun  fix  /,), 


4x  3 

jkZ  nnfi[cxp{—jk3)  . 


^  ■  47r 


sin^i/,. 


in  which 

/«.,  =  /  ^  «cpL;  fc(<y  cos  V*  +  T  cos  fix )]  do' dr' 


or 


...  -  i:{r-Qa-v> 

■exp\jk{oconif}  r  con  fix)]  do' dr' . 


(5.23) 

(5.24) 

(5.25) 


(5.26) 
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The  integrals  J,,,  correspond  to  an  edge  (1)  contribution  to  the  non- 
uniform  field.  The  latter  is  the  field  radiated  by  the  Michaeli’s  fringe  equiv¬ 
alent  current  and  equals 


(5.27) 


with 


4.(co.X-co.^!) 

exp[-jfc(s 


(5.28) 


The  edge  currents  /"“  and  Af"“  have  been  defined  by  eqs.  (5.8)-(5.9)  and 


Figure  5.6:  Geometry  for  the  surface  current  integration. 

(5.10)-(5.12).  The  field  of  eq.  (5.28),  when  superimposed  to  the  Physical 
Optics  contribution  is  equivalent  to  the  GTD  comer  diffraction  recently 
developed  by  Brinkley  and  Marhefka  [2]. 

On  the  other  hand,  the  integrals  if*,,,  corresponding  to  a  tninf.ii.tinn 
effect  of  the  primary  non-uniform  current,  with  the  use  o^  eqs.  (5.2)  and 
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(5.3),  are  written  explicitly  as 


=  (5.29) 


K.  = 


xpf  E  /  i(-i)'y£:;,c«/s!rin{ 

p^i  •'r, 

+flj.  cot  CO.  {]  cot  «■(?)  «'£.  (5.30) 

where  A'(^)  denotes  the  simple  integral 

K{i)  =  [  cxp\jkT'{coB$i-cosfi[)] 

’f  expfjfctr'sin*  +  cos^))  do-'dr',  (5.31) 

which  can  be  trivially  evaluated.  It  can  be  easily  seen  that 

(r(T')  =  r'  sin  a/  nn{0[  —  a)  +  JB,  (5.32) 

where  a  is  the  angle  between  edge  (2)  and  edge  (1)  in  the  counterclockwise 
direction  and  S  is  a  constant.  Hence, 


with 


^^m(0  = 


m)  =  i^,(o  -  MO, 

sin(/3{  -  a)  exp\jk8„  sin*  0[  cos^(^/2)] 


(5.33) 


k’ sin  a  sin* /9(  (^  +  cos  ()(*/  + cos  ^) 

•cxp\jk{i* -f^  -«m  +  «''*rQ^)],  Tn  =  1,2, 


(5.34) 


where  s„,  =  |gm$m|»  A*  bas  been  defined  by  eq.  (5.13),  while 

8in()9J-a)-  . 

1/  =  coeS  =  fi  +  - - •  f -"(cos^i  -  cos^,) 

sin  a  sin  fti 

_  sin  a  sin  0i  cos  4-  coBa(co8^i  —  cob0[) 
sin  a  sin  0[ 


(5.35) 
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Incorporating  eq.  (5.34)  into  eqs.  (5.29)  and  (5.30)  one  obtains 

2  8in(^j  ®)-®*(9n»)  T»rflV/j  1  a!  \ 

"(m)«  - 

•exp[iA!(i‘'-fQ^-s„)l, 


(5.36) 


Kir.).  =  ^  t  W  [yg;(im)W'W(A,  ) 

•cxp(iA:(i‘'-f<,„  -s„)]. 


In  the  above  expressions: 


Sniirj  sin  a 


(5.37) 


0=1 


cot  exp[2jfesn.  cos^(^/2)] 

(/*  +  C08^)(*/ +  cosf)  ’ 

(5.38) 


1  J 

^ — — — E(-i)' 

oniirj  wna  ^ 


/  - 
</r,+r, 


in  ^  cot  exp[2jk9m  co8*({/2)] 


{fl  +  C08()(t/  +  COS() 


(5.39) 


Sriiirj  sin  a 

■tj  - 

^■'ri+n 


cos  ^  cot  cxp[2jk3m  cos^(^/2)] 

(^  +  C08^)(l/  + COS^) 

(5.40) 
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The  pertinent  field  quantities  related  with  the  truncation  effect  read 


jp{2)nu 


and 


j  sinjfi'i  -  g) 
Ink  sin’  0[ 


+ Z  ( 5m )  [cm  A  co«  1  ) 


+  nn  A  cot  !«,*;)]} 

eip[-jt(j  +  Jm  -  ■  fa.)! 


(5.41) 


exp[-jfe(s  +  •  fq^)) 

s 


(5.42) 


It  should  be  noted  that  the  distance  between  Qm  'denoted  by 

s,n-  may  be  *ero  in  case  of  adjacent  edges  or  large  when  and  Qm  '-o  not 
coalesce.  Thus,  the  integrals  nnd  can  be  evaluated  in  closed 

form  when  Sm  =  0  ▼!*  the  Cauchy’s  residue  theorem  and  asymptotically 
for  large  values  of  Sm  via  the  extended  Pauli-Clemmow  method  of  steepest 
descent.  In  the  former  case,  one  obtains 


U{^  -  *') 

sin  a{fi  +  cos  ^)(v  +  cos 

_ sin’/?’ _ 

Til  %\n{0'  —  a)(cos/d  —  cos/?') 

CSC  7  sin  csc^  sin 

cos  ^  -  cos  cos  ^  —  cos 

(5.43) 


H'i’'(8,*^.*') 


U{ir  —  ^')  sin  if>' 
sin  a(/i  4-  cos  ^'){t'  -f  cos  <^) 
sin’  0*  sin 

.f. - L - 21 - 

nj  sin()?'  —  a)(cos/?  —  cos^') 
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.  ]_ 
cos  -  cos  (^)  cos  -  cos  (^)  J  ’ 

(5.44) 

fit  ^•\  -  U{ir -<!>')  COB  4>' 

, _ /9* _ 

Til  8in(/3'  —  a)(co8/9  —  cos/?') 

cot  7  sin  cot  6  sin  (^j^) 

cos  —  cos  cos  ^  —  cos 

(5.45) 

whereas,  a  uniform  asymptotic  expansion  of  these  integrals,  with  an  error 
that  does  not  exceed  0(1//Cm),  with  Km  =  ibs,„Bin*/?,  yields 

«  4  ■  a, 

4ni  \/2irk3„  sin  a  sin  p' 


n>^mc) 


Op»(aw  -  «)(««  -  c(c  -  o„)(c  -  d) 


,  F{^d)  1 

'<KJ-  a„)(<l-c)J’ 

(5.46) 

(5.47) 

(5.48) 

in  which 

c=l  —  /i,  d  =  1  —  Uf 

(5.49) 

and 

Op,  =  1  +  C08[2ni  jVpgVr  -  (-1)^’^'], 

(5.50) 

where  N„ 

integers  most  nearly  satisfying  the  equations 

2n,Ar„ir-(-l)V'  =  (-l)’'-‘»r. 

(5.51) 

One  can,  with  no  difRculty  [40],  obtain  higher  order  terms  of  the  asymptotic 
expansion  of  the  integrals  (5.38)  -  (5.40). 
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5.4  Secondary  non-uniform  edge  equivalent 
currents 

The  analysis  of  the  radiation  of  the  first  order  non-uniform  currents  by 
edge  (1)  merely  included  the  effect  of  truncation  of  these  currents  by  the 
other  edges  of  the  plate.  Of  equivalent  significance  is  the  contribution 
from  a  secondary  effect  which  arises  from  the  non-uniform  currents  that 
are  excited  by  edge  (2)  due  to  the  field  diffracted  by  edge  (1)  that  grazes 
the  plate  and  illuminates  a  section  of  edge  (2),  as  illustrated  in  Fig.  5.7. 


Figure  5.7:  The  secondary  non-uniform  current. 

A  similar  approximation  has  been  developed  by  Sikta  [11],  who  consid¬ 
ered  a  secondary  edge  current  induced  by  the  field  diffracted  by  edge  (1). 
However,  his  development  has  been  based  on  a  modified  version  of  Ryan 
and  Peters  equivalent  currents  [17],  and  though  it  provided  reasonable  re¬ 
sults  for  several  flat  plate  geometries,  it  is  questioned  conceptually  [39]. 
A  natural  modification  of  Sikta’s  idea,  within  the  frame  of  an  asymptotic 
PTD  approach,  is  the  additive  correction  of  the  solution  developed  in  Sec¬ 
tion  3.1  by  the  radiation  of  secondary  non-uniform  equivalent  sources  along 
edge  (2),  which  at  each  point  are  proportional  to  the  tangential  component 
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(to  edge  (2))  magnetic  field  diffracted  by  edge  (1)  and  evaluated  via  UTD. 
Referring  to  Michaeli  [7],  the  non-uniform  current  along  edge  (2)  excited 
»  by  the  UTD  diffracted  field  from  edge  (1)  can  be  written  in  the  form 


/■'■"“(j,)  =  jt;“(ft,ft;ft,Oin,):^i^l^ 

(5.52) 

and  _ 

(5.53) 

where  0;/72j  ^a)  and  m^“(;32,0;/9j,  ^j)  are  determined  from  eqs.  (5.11) 

and  (5.12),  while  the  magnetic  component  H'^  of  the  UTD  diffracted  field 
from  edge  (1)  equals 

-y cxp{-jir/A)  \  \ 

”  =  — — r(  2n. 

$u  (5.54) 

where  [16] 

=  1  -f  cos(2nijV*7r  ±  <f>[) 

(5.55) 

and  integers  most  nearly  satisfying 

2ni  ±  =  tt. 

(5.56) 

It  should  be  also  noted  that  a  factor  of  1  /2  has  been  introduced  in  eqs.  (5.52) 
and  (5.53)  due  to  grazing  incidence. 

Accordingly,  the  field  associated  with  the  secondary  non-uniform  equiv¬ 
alent  currents  defined  by  eqs.  (5.52),  (5.53)  is  given  by 

47rj 

•  [ft  -  ft  cxpljki''  ■  r‘(z,))ifc, 

(5.57) 

where  the  integration  takes  place  along  the  illuminated  (by  the  rays  dif¬ 
fracted  from  edge  (1))  section  of  edge  (2).  The  above  expression  involves 
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Figure  5.8:  Non-parallel  edges  geometry. 


integrals  of  the  form 

/=  /  I  V  H  Hx  vyi./J  exp{jfe{i' .  i\z^)  -  <r(za)  -f-  s**  •  r^(za)]}dza. 

V<r(«a) 

(5.58) 

Although  the  UTD  diffracted  field  from  edge  (1)  is  not  a  valid  approxi¬ 
mation  when  a  corner  is  approached,  because  either  the  edge  itself  is  ap¬ 
proached  or,  at  the  tip,  the  field  structure  is  more  involved  than  a  simple 
superposition  of  singly  and  multiply  diffracted  rays,  the  closed  form  inte¬ 
gration  of  I  is  not  expected  to  yield  a  different  result  from  a  numerical 
integration  along  an  interval  that  excludes  a  small  region  close  to  the  tip, 
due  to  the  regularity  of  the  integrands.  Besides,  closed  form  results  provide 
the  necessary  insight  for  the  assesment  of  the  solution  and  the  possibility 
of  plausible  heuristic  modifications. 

If  the  edges  are  parallel  the  integration  in  eq.  (5.58)  is  trivial.  Let  us 
consider  the  interaction  between  two  non-parallel  edge  as  shown  in  Fig.  5.8 
One  can  integrate  along  the  direction  f  as  shown  in  Fig.  5.6  by  performing 
a  simple  linear  transformation  of  variables.  Since  <r(T)  =  Ar,  where 


A  =  -. 


sin  a 


sin(^;  -  a)’ 
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where 


.  j  — - - - 

sin/Oj  ''  jkAsm^ Pi{a{(f>i)  -  [I  -  I/)] 


|1  —  v\  fi[\\-v\ 


- exp{ j ksm sin^ 0[[a{<l>[)  -  (1  -  i/))}  J ^ — __^exp(-j<*)<fi| 


1  -  1/  <  0 


(5.65) 


One  can  simplify  the  above  expression  in  terms  of  the  edge  transition  func¬ 
tion,  namely 


= 


-j  Bm{P[  -  g) 


”  ifcsinasin/0J  sin^2(l  ~ 

where  =  ksm  sin*  0[  and 


F^iK^,a{<f>[),l-u),  (5.66) 


Oi  —  Oj 


(5.67) 


Substituting  the  result  (5.66)  into  eq.  (5.57)  one  obtains  two  terms  associ¬ 
ated  with  the  end  points  of  the  line  radiation  integral,  namely 


^d,nu  _  ^d,nu  _  j^,nu 


(5,68) 


in  which 


pd,nu  _  pi  („  .  exp(j7r/4)sin(^^  -  a)  sin ^2 
^!^’"‘^2n,(2irib)»/2Bina8in^J(l  -  u) 

•  [4*Wj,^2;/3i,0;n2)  -  ^jmr(/92,^2;)9;,0;n2)] 

•  jcot  i^(«m;a'’^('^i)»l  -  *') 

+  COt  i^(Km;Cr"(^'l).l  -  I')! 

exp[-ife(s,n-s‘'-r,.,)) 


(5.69) 
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When  the  points  Qm  9m  coalesce,  as  in  the  case  of  adjacent  edges 
one  can  repeat  a  similar  analysis  to  identify  a  corner  effect  pertinent  to 
the  secondary  non-uniform  current.  However,  since  the  integrals  involve 
regular  functions  one  can  simply  take  the  limit  of  the  field  in  eq.  (5.69)  as 
Sm  —*  0.  It  can  thus  be  trivially  shown  that  for  Sm  =  0 


i8in(/9J -a)sin^j  . 

4v^;;Ailin 0(1” 

0;  nj)] 

«Cp[-i*(«m  -  i**  •  Tq^)] 


>  ^2;  0;  nj ) 


(5.70) 


where 


■>/^  +  y/i  —  »/’ 


along  with  the  following  selection  of  the  branch: 


(5.71) 


l-y/l  -  t/\  ;  1  -  1/  >  0 

-j\y/l  -  v\  ;1  -  1/  <  0 


(5.72) 


It  should  be  noted  that  the  apparent  singularity  in  eqs.  (5.69)  and  (5.70) 
when  1/  — »  1  is  canceled  by  an  analogous  singularity  of  the  second  term 
composing  the  total  interaction  between  any  pair  of  edges. 


5.5  Numerical  results  and  discussion 

The  effect  of  the  truncation  effect  of  the  primary  non-uniform  current,  out¬ 
lined  in  Sec.  5.3,  as  well  as  the  contribution  of  the  secondary  non-uniform 
current,  presented  iu  Sec.  5.4,  has  been  investigated  in  monostatic  and 
bistatic  cross  section  computations  of  several  plate  geometries  and  they 
have  been  compiled  against  the  comer  diffraction  solution  [2]  as  well  as 
the  Method  of  Moments.  It  should  be  emphasised  that  the  approach  is  not 
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complete  in  many  aspects.  For  instance,  from  a  PTD  point  of  view  particu¬ 
lar  components  of  the  currents  (e.g.,  a  rigorous  tip  excited  current,  an  edge 
wave  current,  currents  excited  by  multiply  bouncing  rays,  etc.)  have  not 
been  included.  Further,  the  predicted  fields  due  to  the  inexact  depiction 
of  the  currents  do  not  satisfy  reciprocity.  However,  the  present  solution 
appears  as  a  computational  efficient  means  to  accurately  approximate  in 
several  cases  the  scattered  fields. 

The  two  levels  of  approximation  of  the  currents  and  the  associated  fields 
are  indicated  as  “asymptotic  solution”  and  “as]rmptotic  solution  +”  at  the 
plots.  In  Figs.  5.9  and  5.10  the  effect  of  the  incorporation  of  the  higher  order 
current  effects  is  evident  away  from  broadside.  It  should  be  noted  that  for 
principal  plane  observation  only  a  ^  polarized  incident  wave  is  expected  to 
excite  higher  order  interactions  between  opposite  edges  of  the  square  plate, 
and,  indeed,  this  is  the  case  where  the  superposition  of  secondary  effects 
improve  the  pattern.  On  the  contrary,  the  comer  diffraction  solution  is 
adequate  in  the  9  polarized  case.  The  minor  discrepancy  when  grazing 
incidence  is  approached  may  be  attributed  to  a  slope  effect  which  has  not 
been  incorporated  at  the  time. 

The  effect  of  the  primary  non-uniform  wedge  current  is  also  clear  in 
bistatic  pattern  calculations  for  the  same  5 A  square  plate.  The  secondary 
current  appears  of  minor  importance  in  this  example,  except  in  the  case  of 
the  cross-section.  When  the  plate  becomes  smaller  and  the  incident 
field  is  closer  to  the  plane  of  the  plate,  the  contribution  from  doubly  dif¬ 
fracted  rays  predicted  via  the  presented  asymptotic  analysis  becomes  more 
significant  as  shown  in  Figs.  5.15-5.18  where  the  bistatic  cross  section  of  a 
2A  square  plate  is  studied. 

Sikta’s  triangle  is  another  example  where  the  method  is  outlined.  Al¬ 
though  the  solution  has  a  remarkable  effect  in  the  sidelobes  close  to  broad¬ 
side  and  when  nose-on  incidence  is  approached,  it  fails  to  predict  the  ap¬ 
preciable  field  at  the  vicinity  of  the  ^  =  150°  angle  for  the  ^  polarization. 
This  discrepancy  may  arise  either  because  of  an  inaccurate  prediction  of 
the  comer  diffracted  field  from  the  back  tip,  cr  due  to  a  strong  tip-to- 
edge  diffracted  ray  (which  is  not  included  due  to  the  two-dimensional  and 
non- reciprocal  nature  of  our  approximation),  or,  finally,  by  travelling  waves 
along  the  edges  of  the  triangle. 

Future  work  on  the  same  level  should  focus  on  a  more  rigorous  develop¬ 
ment  of  the  secondary  non-nniform  current  using  concepts  of  the  Spectral 
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Theory  of  Diffraction.  This  should  enable  the  study  of  slope  effects  and 
should  lead  into  a  more  symmetric  solution.  The  truncation  of  the  sec¬ 
ondary  non-uniform  by  adjacent  corner  is  also  under  investigation.  Further, 
careful  study  of  the  individual  ray  mechanisms  constituting  the  total  solu¬ 
tion  may  enable  heuristic  modifications,  based  on  empirical  arguments,  and 
corrections  in  certain  regions  of  deficiency.  A  more  involved  study  should 
concentrate  in  direct  or  indirect  approaches  towards  the  investigation  of 
the  comer  phenomenon. 
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Figure  5.9:  Principal  plane  RCS  <rtt>  of  a  5A  square  plate. 


Figure  5.10:  Principal  plane  RCS  of  a  5A  square  plate. 


Figure  5.12:  Bistatic  cross  section  of  a  5A  square  plate  for  S'  =  45®  and 
^'  =  0®  at  ^  =  60°. 

95 


Figure  5.17:  Bistatic  cross  section  <t^  of  a  2A  square  plate  for  V 
=  45®  at  ^  =  60®. 


Figure  5.18:  Bistatic  cross  section  of  a  2A  square  plate  for  V  —  60®  and 
=  45®  at  ^  =  60®. 
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Figure  5.19:  RCS  of  Sikta’s  triangle  with  o  =  3A  and  a  =  30°  at 
d  =  90®. 


Figure  5.20:  RCS  <r^/  of  Sikta’s  triangle  with  a  =  3A  and  a  =  30®  at 
e  =  90®. 
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Chapter  6 
Summary 


This  report  summarizes  the  major  theoretical  studies  pertaining  to  the 
scattering  from  flat  plate  structures.  A  new  far  zone  corner  diffraction  co- 
efflcient  has  been  developed  and  tested  against  existing  solutions,  method 
of  moments,  and  measurements.  It  has  been  shown  to  be  very  useful  and 
accurate  for  backscatter  and  especially  for  bistatic  scattering  where  many 
of  the  previous  methods  give  less  accurate  results.  A  new  method  for  de¬ 
termining  an  edge  wave  -  vertex  diffraction  coefficient  for  dipole  excitation 
near  an  edge  has  also  been  developed.  It  is  tested  against  method  of  mo¬ 
ments  and  measurements  with  excellent  results.  This  is  an  intermediate 
step  for  determiiung  a  far  zone  edge  wave  solution.  This  edge  wave  solu¬ 
tion  for  plane  wave  incidence  has  been  approximately  developed  by  relating 
the  first  comer  as  if  it  is  a  source  of  the  edge  waves  and  the  use  of  reci¬ 
procity.  It  is  tested  against  method  of  moments  and  high  resolution  time 
domain  measurements  for  a  plate.  The  higher  order  diffractions  across  the 
face  of  the  plates  has  also  been  developed.  It  includes  comer  to  edge,  edge 
to  comer,  and  comer  to  comer  type  diffraction  mechanisms.  It  is  based  on 
the  truncation  of  the  edge  generated  currents  at  opposite  edges  making  up 
a  plate  stmcture. 

The  theories  in  this  study  have  been  developed  in  such  a  way  as  to  be 
useful  for  efficient  algorithms  to  I  used  in  user  oriented  computer  codes. 
They  have  been  included  in  a  simple  user  oriented  code  called  the  Radar 
Cross  Section  -  Basic  Scattering  Code  (RCS-BSC  Version  2).  It  allows  the 
use  of  multiple  edged  flat  plate  structures  that  can  be  put  together  to  make 
polyhedral  shapes.  It  also  allows  composite  elliptic  cone  frustum  shapes 
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and  ellipsoids  to  model  curved  surfaces.  Monostatic  or  bistatic  results  for 
the  complete  scattering  matrix  are  obtained.  The  input  of  the  shapes  is 
based  on  a  simple  command  word  system. 
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Appendix  A 

Asymptotic  Evaluation  of  an 
Integral  for  Edge  Waves 


We  outline  the  uniform  asymptotic  approximation  of  the  integral 

^-p(*)=/  M-l) 

JO 

for  8?(p)  >  0,  in  terms  of  the  large  parameter  k.  (To  assure  convergence  of 
intermediate  integrals,  let  the  parameter  k  have  a  small  negative  imaginary 
part.)  In  eq.  (A.l),  G{i)  is  a  slowly  varying  regular  function  of  i.  For 
our  purposes,  and  without  loss  of  generality,  it  is  assumed  that  p'(0)  = 
sign(<,)|iy'(0)|,  with  =  0,  and  <  0,  for  every  t.  For  t,  >  0  the 
major  contribution  to  the  integral  arises  from  the  vicinity  of  the  stationary 
phase  point  (t  =  t,)  and  the  end  point,  which  coincides  with  the  branch 
point  singularity  of  the  integrand. 

In  the  particular  case  p  =  1/2,  one  introduces  the  transformation  (  =  y/i 
and  with  H{^)  =  (?(^’),  h{()  —  g{(^)  rewrites  eq.  (A.l)  as 

I-I/M  =  2  r  H{()^xp{jkh{0h  U.2) 

Jo 

where,  now,  h{()  exhibits  three  coUinear,  equidistant  saddle  points  at  (o  = 
0,  (i,2  =  'I'lus  case  is  treated  in  detail  in  [41].  Specifically,  intro¬ 

ducing  the  transformation 


MO  =  oo  -  (a  +  s»)», 


(A.3) 
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with 


<*0  =  M6.a)  = 

a  =  -8ign{<,)|y/i(6.j)  -  M0)l  =  -8ign(<,)|y^(<.)  -^{0)1, 
and  from  the  definition  of  the  parabolic  cylinder  function 

!>-,(*) = — r *’  ^ 

and  the  formula  [42] 

X>_p(a)  =  exp  (jpir)  D-p(-a)  +  exp  {j{p  -  l)ir/2}  Pp_i0z),  {A.5) 

r>-i/a[exp  (iir/4)|o| V^j 

•  -  v^cxp  (-iir/4)  P-i/a[exp  (-yV/4)|a|v^];  <,  >  0 

-4j^^^-i/al«p(i^M)la|V2fcl;  i.  <  0 


(A.6) 


It  should  be  noted  [41]  that,  when  a  is  small, 

<1^1  _  2|ar  2)a|  [  2 

■fa'"”  N  ij-wi  “  \  i«.j»(<.)r  lis"{o)iJ  •  ‘  ’ 

In  the  more  general  case  of  arbitrary  p,  a  rather  heuristic  approach  has 
been  implemented.  For  instance,  the  derivation  of  the  end  point  contribu¬ 
tion  is  obtained  by  introducing  the  expansion 


j(<)«j(0)+j'(0)(  +  2^(*. 


(-4.8) 


The  resulting  canonical  integral  can  be  expressed  in  terms  of  the  parabolic 
cylinder  function  of  order  —p.  Thus, 

^-p{k)  «  f7(0)  r(p)  (2A!)"'/’exp(-y>r/4)  exp{>fc[p(0)  +  oV2]} 

.  M2l!l  ^  I  ®*P  (i»r/4)|o| V^j;  <,  >  0  . . 

.  2|o[  J  \  P_p(exp(iir/4)|o|V^l;  «,  <  0  ‘ 
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